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I.  INTRODUCTION 


There  is  a  continuing  need  for  and  interest  in  simple  models  of  explo¬ 
sively  driven  metal.  Military  applications  have  received  the  most  attention. 


1  2 

beginning  with  the  work  of  Taylor  and  Gurney  on  the  launch  speeds  of  frag¬ 
ments  from  explosively  driven  metal  shells.  Taylor  also  applied  his  ideas  to 

the  formation  of  shaped  charge  Jets^  and  later  published  his  results.1*  Varia¬ 
tions  of  these  basic  ideas  are  still  being  used  today  in  design  and  analysis 

worit.^'^'^  These  ideas  have  also  proven  to  be  useful  in  the  general  field  of 

o 

shock  physics,  and  several  recent  papers  have  appeared  on  applications  of 


G.  I.  Taylor,  "Analysis  of  the  Explosion  of  a  Long  Cylindrical  Bomb  Detonated 
at  One  End,"  paper  written  for  the  Civil  Defence  Research  Committee,  Ministry 
of  Home  Security  in  1941,  available  in  The  Scientific  Papers  of  Sir  Geoffrey 
Ingram  Taylor,  Vol.  Ill,  G.  K.  Batchelor,  ed..  The  University  press, 

Cambridge,  1963*  p.  277. 

2 

R.  W.  Gurney,  "The  Initial  Velocities  of  Fragments  from  Bombs,  Shell, 

Grenades,"  Ballistic  Research  Laboratory  R405,  1943.  (AD  ATI  36210) 

^G.  I.  Taylor,  "A  Formulation  of  Mr.  Tuck's  Conception  of  Munroe  Jets," 
written  in  1943  and  available  in  The  Scientific  Papers  of  Sir  Geoffrey  Ingram 
Taylor,  Vol.  Ill,  ed.  by  G.  K.  Batchelor,  The  University  Press,  Cambridge, 

1963,  P.  358. 

li 

G.  Birkhoff,  D.  P.  MacDougall,  E.  M.  Pugh,  and  G.  I.  Taylor,  "Explosives  with 
Lined  Cavities."  Journal  of  Applied  Physics,  19,  1948,  p.  563* 

^J.  T.  Harrison,  "BASC,  An  Analytical  Code  for  Calculating  shaped  Charge 
Properties,"  Sixth  International  Symposium  on  Ballistics  sponsored  by  the 
Amerioan  Defsnse  Preparedness  Association  (ADPA),  1981,  p.  253* 

^G.  Randers-Pehrson,  R.  R.  Karpp,  C.  E.  Anderson,  Jr.,  and  H.  J.  Blische, 

"SHORTFRAG  Users  Guide,"  Ballistic  Research  Laboratory  Memorandum  Report 
3007,  1980.  (ADB  Qk66kkL) 

7 

L.  R.  Kruse,  "Theoretically  Determined  Fragmentation  Characteristics  and 
Lethality  Estimates  for  a  Newly  Proposed  Tungsten  Alloy  Cased  BLV-97  Combined 

Effects  Bomblet, "  Ballistic  Research  Laboratory  Memorandum  Report  329*1,  1903»(ADC  032497D) 

8 

M.  A.  Meyers  and  L.  E.  Murr,  eds.  Appendix  E,  "Nomograph  for  Determination  of 
Flyer-Plate  Velocity,  11  Shock  Waves  and  Hlgh-Straln-Rate  Phenomena  in  Metals, 
proceedings  of  an  International  Conference  on  Metallurgical  Effects  of 
Hlgh-Straln-Rate  Deformation  and  Fabrication,  Plenum  Press,  NY  and  London, 

1981,  p.  1057. 
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Gurney’s  analysis.  Other  references  to  work  on  this  subject  can  be  found 

in  these  papers.  However,  a  bibliography  compiled  from  these  sources  alone  is 
far  from  complete,  as  will  become  evident  below. 

Here  we  will  express  Taylor's  analysis^  in  simple  formulas  and  extend  it 
to  cases  of  planar  as  well  as  cylindrical  symmetry.  Taylor's  numerical 
analysis  was  presented  in  tabular  form  and  has  been  virtually  ignored.  By 
contrast,  Gurney's  results  were  presented  in  simple  formulas  which  have  been 
widely  used  and  extended.  An  exception  is  Taylor's  simple  formula  for  the 
projection  angle.  This  small  part  of  Taylor's  model  has  also  been  widely  used 
in  a  somewhat  obscure  combination  with  Gurney  formulas  to  obtain  a  projection 
velocity  (magnitude  and  direction).  Here  we  will  develop  time-dependent 
formulas  for  projection  positions  and  velocities  as  simple  as  Gurney  formulas .‘ 
First,  we  will  treat  side-on  cases  in  which  the  detonation  front  propagates 
parallel  to  the  metal  surface  (surface  normal  perpendicular  to  the  propagation 
vector).  Cases  in  which  the  propagation  vector  and  metal  surface  normal  make 
an  angle  between  zero  and  ninety  degrees  can  be  treated  in  the  manner  of 
References  3  and  4.  Finally,  we  will  give  some  simple  formulas  for  head-on 
projection  velocities  (propagation  vector  parallel  to  surface  normal  vector). 

In  Section  II  of  this  report  we  will  establish  the  validity  of  using  the 
exponent  V  =  3  in  the  entropic  equation  of  state.  This  is  required  for  the 
simplification  of  Taylor's  theory  and  its  extension  to  planar  geometry  as 
discussed  in  Section  V.  A  description  and  comparison  of  Taylor's  and  Gurney’s 
models  is  given  in  Sections  III  and  IV  by  way  of  preparation.  Thomas'  model 
is  simplified  and  compared  with  other  models  in  Section  VI.  Comparisons  with 
experimental  data  are  also  made  wherever  possible.  In  Section  VII  we  describe 
Sterne's  extensions  of  the  models  of  his  co-workers,  Gurney  and  Thomas. 

Finally,  we  will  generalize  Sterne's  method  in  order  to  discuss  more  than  one 
layer  of  metal  in  contact  with  more  than  one  layer  of  explosive  in  three 
geometries:  spherical,  cylindrical,  and  planar.  In  particular,  we  will 
develop  simple  design  methods  for  optimizing  certain  types  of  performance. 


q 

G.  E.  Jones,  J.  E.  Kennedy,  and  L.  D.  Bertholf,  "Ballistic  Calculations  of 
R.  W.  Gurney,"  American  Journal  of  Physics,  48.  1980,  p.264. 

10G.  E.  Jones,  "The  Gurney  Equations  for  Multilayered  Fragments,"  Journal  of 
Applied  Physics,  50,  1979,  p.  3746. 

^M.  J.  Kamlet  and  M.  Finger,  "An  Alternative  Method  for  Calculating  Gurney 
Velocities,"  Combustion  and  Flame,  34,  1979,  p.  213. 

12 

E.  W.  LaRocca,  "A  Simplified  Method  of  Calculating  the  Gurney  Constant  of 
Common  Explosives,"  presented  at  the  meeting  of  the  Pyrotechnics  and 
Explosives  Applications  Section  of  the  American  Defense  Preparedness 
Association,  1978. 

13 

E.  W.  LaRocca,  "Advances  in  Predicting  the  Relative  Power  of  High  Explosives 
and  the  Correlation  with  Fragment  Velocity,"  presented  at  the  meeting  of  the 
Pyrotechnics  and  Explosives  Applications  Section  of  the  American  Defense 
Preparedness  Association,  1980. 

14 

D.  R.  Hardesty  and  J.  E.  Kennedy,  "Thermochemical  Estimation  of  Explosive 
Energy  Output,"  Combustion  and  Flame,  28,  1977,  p.  45. 


II.  THE  ADIABATIC  EXPANSION  OF  DETONATION  PRODUCT  GASES 


A  key  element  in  Taylor’s  analysis  was  provided  by  H.  Jones  who  numer¬ 
ically  evaluated  one  of  Taylor’s  integrals.  Jones  worked  with  Taylor  on 
the  Research  Committee  (R.C.)  of  the  Ministry  of  Home  Security  and  wrote 

several  papers  with  him. 15,16, 17  Jones  gave  an  account  of  his  methods  in 

several  Research  Committee  papers  which  are  no  longer  available .18,19,20  por_ 

21  22 

tunately,  his  work  was  later  published.  ' 


In  Reference  21,  Jones  gives  as  his  Equation  (27) 


p  =  Bp 


Y 


(1) 


relating  the  pressure  p  and  density  p  of  the  product  gas  during  its  expansion. 
Here  B  and  y  are  constants.  Equation  (1)  is  called  the  entropic  equation  of 
state  and  applies  to  pol^tropic  gases  for  which  the  internal  energy  depends 

only  on  the  temperature.  For  ideal  gases  at  moderate  pressures  and  tempera¬ 
tures,  y  is  not  much  larger  than  unity.  Jones  estimates  that  y  should  be 
close  to  three  near  the  rear  of  the  reaction  zone  of  a  typical  solid  explosive 
loading  and  remarks  that  such  a  high  value  occurs  because  the  products  ini¬ 
tially  form  a  very  imperfect  gas.  He  derives  the  value  y  =  3  by  assuming  that 
the  ratio  of  the  densities  at  the  front  ( po )  and  rear  (P^)  of  the  reaction 

zone  is  P/p,  s  3/4. 
o  1 


IS 

G.  I.  Taylor  and  H.  Jones,  nNote  on  the  Lateral  Expansion  behind  a 
Detonation  Wave,"  The  Scientific  Papers  of  Sir  Geoffrey  Ingram  Taylor,  Vol. 
Ill,  G.  K.  Batchelor,  ed.,  The  University  Press7  Cambridge, "  1963,  pii  309. 

^G.  I.  Taylor  and  H.  Jones,  "Blast  Impulse  and  Fragment  Velocities  from  Cased 
Charges,"  The  Scientific  Papers  of  Sir  Geoffrey  Ingram ^Taylor,  Vol.  Ill,  G. 
K.  Batchelor,  ed.,  The  University  Press,  Cambridge,  1963,  p.  363. 

17 

G.  I.  Taylor  and  H.  Jones,  "The  Bursting  of  Cylindrical  Cased  Charges,"  The 
Scientific  Papers  of  Sir  Geoffrey  Ingram  Taylor,  Vol.  Ill,  G.  K.  Batchelor, 
ed..  The  University  Press,  Cambridge,  1963,  p. 379. 

1 8 

H.  Jones,  Research  Committee  Paper  166,  no  longer  available. 

19  ' 

H.  Jones,  Research  Committee  Paper  212,  no  longer  available. 

20 

H.  Jones  and  A.  R.  Miller,  Research  Committee  Paper  306,  no  longer 
available. 

21 

H.  Jones,  "A  Theory  of  the  Dependence  of  the  Rate  of  Detonation  of  Solid 
Explosives  on  the  Diameter  of  the  Charge,”  Proceedings  of  the  Royal  Society 
of  London,  189A,  1946,  p.  415.  —  — —  —  ~~  “ 

22 

H.  Jones  and  A.  R.  Miller,  "The  Detonation  of  Solid  Explosives:  the 
Equilibrium  Conditions  in  the  Detonation  Wave-Front  and  the  Adiabatic 
Expansion  of  the  Products  of  Detonation,”  Proceedings  of  the  Royal  Society 
of  London,  194a,  1948,  p.  480. 

23 

R.  Courant  and  K.  0.  Friedrichs,  Supersonic  Flow  and  Shock  Waves, 
Interscience,  NY,  1948,  p.  6. 
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In  Reference  22  Jones  and  Miller  used  the  following  virlal  equation  of 
state 

2  3 

pV'/N  =  RT  +  bp  +  cp  +  dp  (2) 

to  calculate  pressure  and  volume  (V*)  of  N  moles  of  gas  at  temperature  T  re¬ 
sulting  from  the  detonation  of  one  mole  of  explosive.  They  also  calculated 
product  species  for  TNT  and  predicted  an  approximately  linear  dependence  of 
both  detonation  velocity  and  detonation  pressure  on  loading  density.  In  Part 

•2 

B  of  their  paper,  they  give  numerical  examples  for  TNT  at  a  high  (1.5  g/cm  ) 

"2  W 

loading  density  and  a  low  (1.0  g/cnr)  loading  density.  For  the  high  density 
loading,  their  Table  III  gives  the  pressure,  temperature,  inverse  density,  and 
effective  Y  value  calculated  for  the  entire  expansion  process.  Here  we  are 
particularly  Interested  in  the  y  values.  At  the  detonation  pressure  achieved 

at  the  rear  of  the  reaction  zone  (p^  =  15.88  x  1010  dyne/cm2)  they  find  y  = 

3-36  (and  P  =  1.95  g/cm^).  By  the  time  the  pressure  has  dropped  to  6.25  x 

10  dyne/cm  ,  p  =  0.65  g/cmJ  and  Y  =  2.39.  Between  these  two  points  we  note 
that  a  y  value  obtained  by  taking  half  of  (3-36  +  2.39)  is  2.9.  Their  Table 

continues  down  to  p  =  2.818  x  10^  dyne/cm2  and  P  =  0.003  g/cm^  with  7  =  1.27. 
Their  Table  V  gives  similar  results  for  the  lower  loading  density  with  y  rang¬ 
ing  from  2.113  at  the  detonation  pressure  down  to  1.20  near  atmospheric  condi¬ 
tions.  Normally  cast  TNT  for  military  applications  ha3  a  loading  density 

•3 

slightly  greater  than  1.6  g/cra  so  that  a  suitable  effective  value  over  the 

range  of  Interest  to  metal  acceleration  (2.1  >  P  >  0.6  g/cmJ)  ought  to  be  y  = 
3*  Jones  and  Miller  go  on  to  compare  their  calculated  species  concentrations 
favorably  with  experiment,  Implying  that  their  p,  T,  p  and  7  estimates  ought 
to  be  approximately  correct  also,  even  though  no  direct  comparison  of  these 
quantities  with  experimental  measurements  Is  possible  even  today.  These 

2H 

results  of  Jones  and  Miller  were  again  used  by  Taylor  In  a  later  paper  in 
which  he  demonstrated  for  the  first  time  the  possibility  of  a  spherical 

detonation  wave,  contradicting  an  opinion  previously  expressed  by  Jouget  .25 

Many  other  authors  have  followed  Jones  by  using  y  =  3  to  describe  the  high 
pressure  gases  which  result  from  the  detonation  of  condensed  explosives  near 

pr 

their  maximum  loading  density.  Jacobs^  has  pointed  out  that  both  the 

27 

Klstiakowsky-Wllson  equation  of  state  and  that  used  by  Jones  and  Miller 

2^G.  I.  Taylor,  "The  Dynamics  of  the  Combustion  Products  behind  Plane  and 

Spherical  Detonation  Fronts  in  Explosives, "  Proceedings  of  the  Royal  Society 
of  London,  200A,  1950,  p.  235. 

2^M.  Jouget,  C.  R.  Academy  of  Science,  1M,  Paris,  1907,  P-  633* 

26  r  — 

S.  J.  Jacobs,  "Recent  Advances  in  Condensed  Media  Detonations,"  American 

Rocket  Society  Journal, ‘30.  I960,  p.  151. 

27 

G.  B.  Klstiakowsky  and  E.  B.  Wilson,  Jr.,  OSRD  Report  1lH,  US  National 
Defense  and  Research  Committee  of  the  Office  of  Scientific  Research  and 
Development,  19JJI. 


10 


(Equation  (2)  above)  lead  to  calculated  results  which  can  be  expressed  by 
Equation  (1)  above  with  7=3  down  to  one  per  cent  of  the  detonation  pressure, 
with  an  accuracy  acceptable  for  hydrodynamic  usage. 

In  1968,  Jacobs,  together  with  Kamlet,  proposed  a  simple  method  of  calcu¬ 
lating  detonation  properties,  that  is,  properties  at  the  rear  of  the  reaction 

pO 

zone.  This  was  the  first  in  a  long  series  of  papers  which  has  not  yet  been 
29-35  33 

completed.  The  sixth  paper  in  this  series  gives  an  empirical  expres¬ 

sion  for  7  in  terras  of  the  loading  density  with  7  again  approximately  equal  to 
three  for  military  explosives  at  typical  loading  densities.  In  connection 

with  this  paper,  we  also  note  a  recent  paper  by  Andersen 

32 

The  fifth  paperJ  in  the  above  series  provides  us  with  a  table  of  pres¬ 
sures  and  inverse  densities  (specific  volumes)  for  various  conditions  during 
the  expansion  of  the  products  from  the  detonation  of  twelve  common  explosive 
loadings,  calculated  by  using  a  Kistiakowsky -Wilson  equation  of  state: 


pV7(NRT)  »  1  +  Xe*xf 


(3) 


28 

M.  J.  Kamlet  and  S.  J.  Jacobs,  "Chemistry  of  Detonations,  I:  A  Simple 
Method  for  Calculating  Detonation  Properties  of  C-H-N-0  Explosives ,"  journal 
of  Chemical  Physics,  48,  1968,  p.  23. 

29 

M.  J.  Kamlet  and  J.  E.  Ablard,  "Chemistry  of  Detonations,  II:  Buffered 
Equilibria,"  Journal  of  Chemical  Physics,  48,  1968,  p.  36. 

30 

M.  J.  Kamlet  and  C.  Dickinson,  "Chemistry  of  Detonations,  III:  Evaluation 
of  the  Simplified  Calculational  Method  for  Chapman-Jouget  Detonation 
Pressures  on  the  Basis  of  the  Available  Experimental  Information,"  Journal 
of  Chemical  Physics,  48,  1968,  p.  43. 

3 1 

M.  J.  Kamlet  and  H.  Hurwitz,  "Chemistry  of  Detonations,  IV:  Evaluation  of  a 
•Simple  Predictional  Method  for  Detonation  Velocities  of  C-H-N-0  Explosives," 
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a  t 

where  X  =  KZx-^k^/l  V'  (T+G)  3  with  x^  the  mole  fraction,  kA  the  covolume  fac¬ 

tor,  and  a',  (8,  K  and  $  empirical  constants.  If  we  use  this  table  to  plot  p 

3 

versus  p  ,  we  see  that  Equation  (1)  with  7=3  can  be  fitted  very  closely  to 
these  calculated  results,  at  least  over  the  range  of  p  of  Interest  to  metal 

26 

acceleration,  confirming  Jacobs’  remark.  Figure  1  illustrates  this  for  two 
explosive  loadings  which  are  not  very  different  from  ones  used  in  military  ap¬ 
plications.  The  values  marked  X  and  0  were  calculated  from  their  table. 

3 

Kamlet  and  Hurwitz  also  quote  values  for  TNT  at  P0  =  1.55  g/cm  which  are 

3  22 

quite  close  to  those  calculated  for  =  1.50  g/cm  by  Jones  and  Hiller, 
using  a  very  different  equation  of  state  (Equation  (2)  versus  Equation  (3)). 

We  conclude  that  Equation  (1)  with  7  =  3  is  a  satisfactory  approximation  for 
most  applications.  In  fact  it  is  much  better  than  we  might  expect  from  the 
estimates  given  by  Jones  and  Miller,  quoted  above.  At  present  we  have  no 
satisfactory  explanation  for  this  remarkable  observation.  The  fact  that  the 
range  .6  <  p  <  2  is  centered  near  p  =  1  where  p  is  independent  of  7  may  be  of 

some  help,  but  this  does  not  explain  the  agreement  over  the  range  0.2  <  P^  < 

8. 


If  we  adjusted  7  as  well  as  B  in  Equation  (1)  we  might  obtain  slightly 
better  agreement.  However,  there  is  little  justification  for  doing  this  since 
the  values  being  fitted  can  be  calculated  from  various  equations  of  state 
chosen  for  convenience  rather  than  theoretical  reasons.  Since  almost  any 
reasonable  form  will  do,  we  should  not  place  much  weight  on  the  fact  that  a 
particular  form  gives  us  reasonable  results.  As  Jones  and  Hiller  remarked 
22 

long  ago,  we  have  little  information  to  guide  us  in  choosing  a  correct  form 
for  the  equation  of  state.  There  is  no  virtue  in  using  complicated  empirical 
equations  which  are  designed  to  cover  critical  phenomena  controlled  by  weak 
molecular  attractions.  The  best  that  can  be  done  is  to  choose  the  simplest 
equation  which  is  adequate  for  one’s  purpose,  as  they  did  and  as  we  shall  do. 

If  we  let  7  =  3  in  Equation  (1),  analytical  simplifications  become  possi¬ 
ble,  as  we  shall  see  presently.  One  example  of  such  simplification  has  already 

37 

been  given  by  Aziz  and  co-workers.  who  compared  the  results  of  an  exact  so¬ 
lution  for  rigid  piston  loading  (possible  if  7  =  3)  with  numerical  solutions 
in  the  parameter  range  2.5  <_  7  <  3.5.  They  concluded  that  only  a  one  percent 
error  would  be  made  in  calculating  the  energy  transmitted  to  the  piston,  in 
spite  of  variations  in  the  detonation  pressure  and  velocity  of  about  fifty  per 
cent.  Similar  conclusions  were  reached  earlier  by  Gurney  in  his  discussion  of 
78 

gas  leakage. 


37 
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— 

R.  W.  Gurney,  "Fragmentation  of  Bombs,  Shells  and  Grenades,"  Ballistic 
Research  Laboratory  Report  635,  1947.  (ADB  800451) 
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Figure  1.  Linear  Relation  between  Pressure  and  Density  Cubed. 


III.  TAYLOR'S  ANALYSIS  OF  A  TUBULAR  BOMB 


Since  the  details  of  this  model  are  not  widely  known,  we  will  summarize 
them  here.  Figure  2  shows  a  longitudinal  cross  section  of  an  infinitely  long 
tubular  bomb  in  the  neighborhood  of  a  planar  detonation  front  (solid  vertical 
line)  which  is  moving  to  the  right  at  constant  speed,  D.  To  an  observer  in 
the  laboratory,  the  process  of  gas  expansion  and  metal  acceleration  appears  to 
be  complicated,  since  the  gas  first  moves  forward  following  the  detonation 
wave  front  and  later  moves  backward  away  from  it  so  that  it  acts  twice  on  the 
same  element  of  the  case,  as  we  shall  see.  However,  to  an  observer  who  trav¬ 
els  with  the  detonation  front  and  looks  back  at  the  expansion  process,  these 
complexities  disappear,  since  the  gas  pressure  and  velocity  at  a  given  dis¬ 
tance  behind  the  front  (dashed  vertical  line)  do  not  change  with  time  (steady 
state  condition).  Moreover,  the  velocity  of  a  case  element  along  its  length 
ds  is  constant  and  equal  to  D,  since  there  is  no  stretching  along  ds.  The 
transformation  between  the  laboratory  coordinate  X  and  the  coordinate  x  used 
by  the  moving  observer  is 


X  =  -  (x-Dt)  +  XQ,  (4) 

as  shown  in  Figure  2.  Here  t  is  the  time,  which  is  taken  to  be  the  same  for 
both  observers,  since  detonation  wave  speeds  are  negligible  compared  to  the 
speed  of  light.  The  radial  coordinates  R  =  r  are  also  the  same  for  both  ob¬ 
servers  as  is  the  (lgnorable)  equatorial  angle  measured  about  the  x-axis.  The 

time  derivative  of  Equation  (4)  gives  us  the  relation  between  the  speeds  U  =  X 

and  u  =  x,  using  the  dot  convention  for  time  derivative, 

U  =  D  -  u.  (5) 


Let  M  be  the  constant  mass  per  unit  length  of  the  cylindrical  shell  casing 
and  r  be  the  radial  position  of  the  center  of  a  ringlike  element  of  the  case 
wall,  namely,  the  average  of  the  inner  and  outer  radii  of  the  casing  at  any 
time  t  with  r  much  larger  than  the  wall  thickness.  If  we  take  the  ring  to  be 
of  unit  length,  then  the  area  acted  on  by  the  gas  pressure  is  2wr  times  unity 
and  the  mass  per  unit  area  is  M/(2irr).  The  pressure  always  acts  normal  to  the 

2 

metal  surface  and  balances  the  centrifugal  force  [M/(2xr)]  (D  /R  ),  where  R 

c  c 

is  the  radius  of  curvature  (see  Figure  2),  the  inverse  of  the  curvature  of  the 
arc  ds.  Thus, 


M_ 
2t rr 


(6) 


where  we  have  retained  the  exact  expression  for  the  curvature.  For  a  suffi¬ 
ciently  heavy  case,  tan2  $  = /^l\2  will  be  much  less  than  unity  (see  Figure 

Vdxf  o 

2),  approaching  unity  only  as'  4>  approaches  45  .  Thus  it  is  a  good  approxima¬ 
tion  to  neglect  it  in  the  denominator  of  Equation  (6)  as  Taylor  implicitly 
did.  However,  there  is  no  difficulty  introduced  by  retaining  it,  at  least  for 
a  while,  so  we  will. 
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Figure  2.  Tubular  Case  Expansion  Geometry. 


If  we  multiply  both  sides  of  Equation  (6)  by  (MflrdrJ/MD^)  =  (4irr  tan  dx) 
/(MD2),  we  find 


tan  <f>  d(tan  cf) j 
[l  *  tan2  J3/2 


2  sin  cf>  d  pd(r2) 

MD2 


(7) 


2  2  2 

where  we  have  used  the  identities  d(tand>)  =  sec  4>d<p,  1  +  tan  4>  3  sec  and 

2 

2  rdr  =  d(r  ).  The  left  side  of  Equation  (7)  is  readily  integrated  from <p  =  0 
to  obtain 

2(1  -  cos  <f>}  =  4  sin2(<J>/2)  (e) 


by  another  Identity.  If  we  approximate  Equation  (8)  by  retaining  only  the 
first  order  term  in  a  binomial  expansion. 


2(1  -  cos  <J>) 


(l  -  -  1--r 

)  =  2  [l  -  (1  +  tan2  pyh\ 

l  sec  p 

J  L 

(9) 

\ 

r. 

'  2  .  .-i  2. 

we  recover  Taylor's  approximation.  The  same  thing  can  be  accomplished  by  re¬ 
taining  only  first  order  terms  in  series  expressions  of  either  side  of  Equa¬ 
tion  (8).  Thus  2(1-  cos<j>]  ~  2  [  1  -(1-^2)]  =  <t>2  =s  tan20,  or  4  sin2(d>/2)  ss 
4(4>/2)2  =  <t>2  ss  tan2d>  for  small  <p . 


We  are  interested  in  the  square  root  of  Equation  (8).  In  moving  coordi¬ 
nates,  the  components  of  the  case  velocity  are  vp  =  dr/dt  =  D  sin  4>  and  v^  = 

dx/dt  =  D  cos  <f>  from  Figure  2.  In  laboratory  coordinates  the  components  are 
V  =  v  and  V  =  D(1  -  cos  <t> ) ,  so  the  magnitude  of  the  case  velocity  in 

.  ^  I 

laboratory  coordinates  is 


=  Sin2  cf>  +  (1  -  cos  p)2  =  D^2(l  -  cos  <J>)  ( 

=  2D  sin  (p/2)  , 


10) 


2  2 

where  we  have  used  the  identity  sin  <f>  +  cos  4>  =  1  and  Equation  (8).  Equation 
(10)  is  the  only  formula  in  Taylor's  analysis  which  has  received  wide  use, 
clearly  because  of  Its  simplicity.  We  also  note  by  another  identity  that 


V  1  -  COS  P 

T7—  “  - - - - —  =  tan  (<p/ 2) 

V  sm  d>  Y/ 

r  Y 


(11) 


so  that  the  angle  of  projection  in  laboratory  coordinates  is  p/2.  Equations 
(10)  and  (11)  are  exact.  If  we  follow  Taylor  and  use  tan  </>  cs  2  sin(4>/2)  for 
the  integral  of  the  left  side  of  Equation  (7),  we  make  less  than  a  5$  error 

for  <  20°.  The  exact  expression  in  Equation  (8)  has  already  been  used  by 
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39  40 

others.  In  the  journal  literature  It  was  used  by  Allison  and  co-workers  ’ 
who  wrote  dtf>/ds  =  (dr/ds)  (d<£/dr)  =  sin4>(d$/dr)  for  the  curvature  in  Equation 
(7).  They  did  not  point  out  that  this  is  equivalent  to  retaining  the  exact 
curvature  instead  of  approximating  the  curvature. 

In  order  to  integrate  the  right  side  of  Equation  (7)  let  us  write  down  a 
few  more  useful  relations.  For  an  explosive  gas  mixture  in  a  rigid  cylinder, 
such  as  a  hydrogen/oxygen  mixture  in  a  shock  tube  experiment,  the  equation  of 
mass  conservation  in  moving  coordinates  is 


dm  =  pQAQDdt  =  pAudt 


(12) 


where  Ddt  and  udt  are  the  original  and  compressed  length  of  a  gas  volume  ele¬ 
ment  of  mass  dm  and  constant  cross  sectional  area  A  =  A  .  Here  p  is  the  un- 

o  o 

compressed  gas  density,  while  P  is  its  compressed  density.  Since  the  tube  is 
rigid,  A  is  constant  and  PqD  =  pu.  However,  for  the  gaseous  detonation  prod¬ 
ucts  of  a  solid  explosive  high  density  loading  expanding  in  a  metal  tube,  A 
will  not  remain  constant.  For  a  tube  with  a  circular  cross  section,  the  area 

2  2 

before  expansion  is  AQ  =  w rQ  ,  while  after  some  expansion  it  becomes  A  =  * r  > 

A  .  If  we  use  these  expressions  in  Equation  (12)  and  cancel  we  find 
o 


PQDrQ^  =  pur^  (13) 

for  mass  conservation.  Now  pQ  is  the  loading  density  of  the  solid  explosive 

before  detonation.  Since  steady  state  flow  conditions  hold  in  moving  coordi¬ 
nates,  the  equation  of  continuity  is 


i  3?  (rpv)  *  97  tpuJ  *  0 


(14) 


where  v  is  the  radial  component  of  the  gas  velocity  (as  contrasted  with  v^, 

the  radial  component  of  the  metal  velocity).  In  Table  1  below  we  will  note 
that  v  <<  u  «  D,  so  that  the  first  term  in  Equation  (14)  might  be  neglected  on 
these  grounds  for  r  values  of  Interest.  This  would  make  (pu)  independent  of  x 
as  it  is  in  Equation  (13).  However,  it  is  more  Interesting  to  conclude  from 
Equation  (13)  that  (pu)  is  independent  of  x,  so  that  (rpv)  in  Equation  04) 

2 

must  be  independent  of  r.  Since  P  is  inversely  proportional  to  r  from 
Equation  (13)  when  u  «  D,  then  v  must  be  proportional  to  r  under  these  condi¬ 
tions.  Thus  we  may  write 


V  *  r(VB/rB)’ 


(15) 


39 

F .  E.  Allison  and  j.  T.  Schriempf,  "Explosively  Loaded  Metallic  Cylinders, 
II,"  Journal  of  Applied  Physics,  31#  I960,  p.  846. 

40 

F.  E.  Allison  and  R.  W.  Watson,  "Explosively  Loaded  Metallic  Cylinders,  I," 
Journal  of  Applied  Physics,  31#  I960,  p.  842. 
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Table  1.  Taylor’s  Results  for  TNT 
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where  v  and  r  are  constants  which  we  will  define  later.  Equation  (15)  was 
B  B  ^ 

Introduced  by  Gurney^  as  a  postulate,  although  it  follows  from  the  equation  of 
continuity  under  the  conditions  mentioned.  Gurney  gave  no  Indication  that  he 
ever  noticed  this  fact.  Instead,  his  assumption  had  quite  a  different  motiva¬ 
tion  as  we  shall  see. 

In  the  Chapman-Jouget  model  of  detonation,  the  approximation  is  made  that 
all  chemical  reaction  is  completed  in  such  a  short  time  that  the  width  of  the 
reaction  zone  is  negligible.  Thus  we  may  approximate  changes  in  variables 
like  p,  p  and  u  by  discontinuous  jumps  Instead  of  using  derivatives  with  very 
large  values.  If  we  use  the  subscript  unity  to  denote  values  just  behind  this 
very  narrow  reaction  zone  and  the  approximation  r1  =  rQ  (no  expansion  in  the 
reaction  zone),  then  Equation  (13)  becomes  1 

poU  =  P1U1 

as  in  a  rigid  tube.  Similarly,  we  can  use  Newton's  second  law  in  its  impulse- 
momentum  form  to  derive  another  jump  condition.  The  impulse  delivered  across 
the  reaction  zone  is  (p^  -  pQ)Adt  which  is  equal  to  the  change  in  momentum 

across  this  zone,  namely, dm(D  -  u^  where  dm  is  given  by  Equation  (12).  If  we 

neglect  the  (ambient)  pressure  p  at  the  front  of  the  reaction  zone  compared 

o 

to  the  detonation  pressure  p^  at  its  rear  and  cancel  Adt,  we  obtain 

P1  =  pqD(D  -  u^*  (17) 

We  need  one  other  relation  in  order  to  integrate  the  right  side  of  Equa- 

4  1 

tlon  (7).  This  ia  the  strong  form  of  Bernoulli's  law  for  steady  flow 

f9\  (“2  *  y2)  ~  15  lUl2  ’  “2)  (18) 

2  2  2 

since  v  <<  u  as  D  as  mentioned  above.  The  approximate  differential  form  of 
Equation  (18)  is  then 


dp  :  -  pudu.  ( 19 ) 

Every  science  student  is  familiar  with  the  integral  of  Equation  (19)  in  the 
case  of  an  incompressible  fluid  for  which  p  =  p  for  any  p  and  u  (in  the  ab- 

°  2  2 

sence  of  external  fields  like  gravity),  namely,  p  +  ^u  =  p^  +?2^ui  *  *n 
the  present  case,  of  course,  p1  >  pq. 

.  2 

The  variable  on  the  right  side  of  Equation  (7)  was  chosen  to  be  r  so  that 
we  may  use  Equation  ( 1 3 )  to  introduce  the  new  variable  1/(pu).  Then  the  right 
side  of  Equation  (7 )  becomes 


4  1 


R.  Courant  and  K.  0.  Friedrichs,  Supersonic  Flow  and  Shock  Waves, 
Interscience,  NY,  1948,  p.  22. 
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2up  r 
o  o 
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MD 


=  f-E.  +  u  -  d1 
aD  Lpu  J 


=  4  sin  (<J>/2) 


(V/D)2J»  tan2  (J > 


(20) 


1 


which  is  Taylor's  Equation  (11)  when  tan  <j>  Is  used.  In  the  second  line  of 
Equation  (20), we  have  Introduced  Taylor's  notation  for  the  metal  case  mass  to 
explosive  charge  mass  ratio 


a 


M 


(p  TTr  ) 
o  o 


M 

C 


(21) 


where  C  Is  the  explosive  mass  per  unit  length  of  cylinder.  We  have  also  Inte¬ 
grated  by  parts  from  the  highly  compressed  state  1  to  an  expanded  state  some 
distance  x  behind  the  detonation  front  (zero  width  reaction  zone).  In  the 
third  line  of  Equation  (20),  we  have  used  Equation  (19)  to  evaluate  the 
remaining  Integral.  In  the  fourth  line,  we  have  eliminated  p  D  in  Equation 
(17)  by  using  Equation  (16)  to  obtain  ° 


plul 


+  u1  =  D 


(22) 


then  used  Equation  (22)  In  Equation  (20).  In  the  last  two  lines  we  have  used 
Equations  (8),  (9),  and  (10)  to  complete  the  derivation.  Equation  (20)  gives 
V  as  a  function  of  the  variables  p,  p  and  u  and  the  parameters  D  and  a  . 

Taylor  went  on  to  find  the  shape  of  the  case  by  numerically  Integrating  the 
equation 


x 


dr/ tan  $ 


which  is  derived  from  tan<£=  dr/dx. 


(23) 


Equation  (20)  is  rarely  cited  in  the  literature  since  p  and  u 
pear  as  simple  functions  of  r,  although  (pu)  has  a  simple  form  in 
(13).  Instead,  u  must  be  found  from  Equation  (18)  with  u^,  given 

(16).  As  was  mentioned  above,  H.  Jones  evaluated  the  integral  in 


do  not  ap- 
Equatlon 
by  Equation 

Equation 
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3  5 

(18)  for  Taylor.  Once  this  was  done  (using  =  1.51  g/cm  and  D  =  6.38  x  10 
cm/s),  Taylor  presented  his  solution  in  tabular  form.  Table  1  here  is  an 
abbreviated  version  of  his  results  with  one  row  added  (for  u  =  D)  and  some 
additional  columns — two  giving  tf>  in  degrees  for  the  exact  and  approximate 
expressions  in  Equation  (20),  using  the  example  «  =  2.  From  this  table  we  see 
that  by  the  time  l  =  r/rQ  =  1.5,  V  has  already  achieved  90$  of  the  value  it 

would  achieve  if  the  case  were  ductile  enough  to  avoid  fragmentation  until 

r/r  =2.  Calculations  for  large  values  of  r/r  are  purely  academic,  since  no 
o  o 

casing  material  will  stretch  indefinitely.  This  is  emphasized  by  the  fact 

that  the  angle  <f>  found  by  using  the  exact  curvature  becomes  imaginary  as  soon 

as y[W  (V/D)  >2.  This  corresponds  to  the  physical  fact  that  V  cannot  exceed 

D.  The  value  <f>  =  80.97°,  found  by  using  the  approximation  tan2^,  is  without 

significance  since  the  approximation  ceased  to  be  valid  long  before 

r/r  -r™3  .  Since  practical  metal  cases  usually  burst  near  r/r  =  1.5,  we  see 
o  o 

that  the  range  of  pressure  and  density  which  can  be  adequately  represented  by 

Equation  (1)  with  y  =  3  is  covered  quite  well.  As  can  be  seen,  approximate 

values  are  very  close  to  exact  values  in  this  range.  Finally,  we  note  that  u 

is  Initially  somewhat  smaller  than  D,  then  becomes  slightly  larger  than  D  for 

r/rQ  >  1.2.  If  we  recall  Equation  (5),  it  becomes  clear  that  to  a  laboratory 

observer  the  gas  tirst  appears  to  move  forward  after  the  detonation  front  (U  > 

2 

0),  then  backward  (U  <  0)  for  any  given  X  position.  The  approximation  v  << 

2  2 

u  ss  D  also  becomes  clear,  since  v  <  Vs s  0.2D  for  reasonable  a  «  2  in  the 
fifth  column.  ” 

Taylor  did  not -have  cased  TNT  observations  to  compare  with  his  calcula¬ 
tions  when  he  wrote  his  paper,  so  he  compared  his  V  values  with  observations 
made  for  tetryl  explosive  cased  in  steel  and  found  rough  agreement.  Tetryl  is 
similar  to  TNT  since  it  too  has  an  hexagonal  trinltro  ring.  However,  it  has 
an  extra  nitro  group  bonded  to  its  methyl  group  via  a  nitrogen  atom,  so  we 

1 6  17 

might  expect  some  disagreement.  Taylor  later  collaborated  with  H.  Jones  ' 
and  compared  calculated  case  shapes  and  velocities  with  experimental  values 
for  composition  B  explosive  (60/40  RDX/TNT)  encased  in  steel,  using  Jones’ 
calculations  for  this  mixture.  They  found  agreement  over  the  range  observed, 

namely,  0*6  <  o  <  13.5,  which  corresponds  to  26.5°  >  <f>  >  6.5°  near  bursting. 

In  recent  years  cylinder  tests  have  become  standardized  as  explained  for 

34 

example  by  Short,  who  gives  measured  cylinder  wall  velocities  as  a  function 

■3 

of  case  expansion  for  two  densities  of  TNT  ( p  =  1.45  g/cm  and  P  =  1.63 
^  00 

g/cnr),  which  bracket  the  densities  being  used  by  Taylor  in  his  calculations 

(P  =  1.51  g/cm  ).  For  the  standard  copper  cylinder, a  =  M/C  =  (pm  /  p  ) 
o  2  3  i*»  o 

^rout/rin^  *  1]  =  4.026/ PQ,  since  p Cu  =  8.92  g/cm  ,  the  outer  cylinder 

radius  r  .  =  1.53  cm  =  r  and  the  inner  radius  r,  =  1.27  cm.  For  p  =  1.51 
t  out  o  in  o 

g/cnr  we  find  a  =  2.67.  We  can  use  this  value  together  with  D  =  .6.38  mm/^ts 

to  find  V  (mm/M^  =  5.52  times  the  values  in  column  5  of  Table  1.  These  values 
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are  about  1 0%  lower  than  the  experimental  values  in  column  10,  consistent  with 
Taylor’s  neglect  of  the  radial  gas  velocity  component,  an  approximation  which 
becomes  worse  as  £  increases.  The  experimental  values  were  obtained  by 

plotting  Short's  values^  versus  r/rQ  =  1  +  (r  -  r  )/rQ,  since  his  expansion 

variable  is  (r  -  r  ),  followed  by  a  linear  interpolation  between  the  cases  he 
0  1  3 

gives  for  P  =  1.45  g/cm  and  P  =1.63  g/cm  .  The  method  of  interpolation  is 
o  o 

unimportant,  since  the  values  in  column  9  lie  below  the  lower  curve  for  p  = 

3  ° 

1.45  g/cm  instead  of  between  the  two  curves.  We  can  imagine  a  variety  of 

corrections  for  both  the  calculated  and  measured  values.  For  example,  we 

might  use  a  more  modern  value  of  D  closer  to  7  mm/MS,  as  indicated  in  the 

survey  by  Kamlet  and  Hurwitz.^  We  might  also  use  the  higher  pressure  values 
•  22 

given  by  Jones  and  Miller  and  adjust  u  accordingly.  On  the  experimental 
side  we  might  correct  for  the  neglect  of  Taylor’s  angle  in  Short's  values 
(which  are  V  ,  not  V)  by  eliminating  4>  between  Equations  (10)  and  (11)  to 
obtain 

V  =  D  £ 2  jl  -Vl  -  (Vr/D)2J]  ^  (24) 

However,  this  leads  to  corrections  only  in  the  third  decimal  place.  Even 
corrections  for  observing  the  outside  instead  of  the  center  of  the  thinning 
case  wall  are  beyond  the  accuracy  of  the  experiment,  since  the  data  reduction 
process  includes  numerical  differentiation  of  displacement  versus  time  curves. 
Altogether,  the  agreement  is  remarkable. 


IV.  GURNEY’S  MODEL  AND  THOMAS'  SYNTHESIS 

2 

It  is  clear  from  Section  3  of  his  report  that  Gurney  wanted  to  explain  an 
experimental  observation;  namely,  the  fact  that  the  fragment  launch  speeds  of 
very  different  size  weapons  containing  the  same  explosive  seemed  to  depend 
mainly  on  the  ratio  of  the  mass  of  the  explosive  to  the  mass  of  the  metal; 

that  is,  (C/M)y,  where  y  is  a  power  near  0.22  for  large  bombs,  but  closer  to 
0.50  for  small  projectiles.  This  led  him  to  his  basic  assumption  which  he 
stated  in  his  abstract  and  repeated  in  Section  4  of  his  report:  namely,  the 
contribution  made  to  the  total  kinetic  energy  by  the  detonation  of  each  unit 
mass  of  explosive  is  independent  of  the  size  of  the  projectile.  He  used  the 
symbol  E  to  denote  this  constant  energy  contribution  per  unit  mass  and  C  =  pQ 

rQ2  for  the  explosive  mass  per  unit  cylinder  length.  Thus, EC  is  the  energy 

contribution  per  unit  length.  Initially,  the  energy  released  by  the  detona¬ 
tion  appears  entirely  as  the  internal  energy  of  the  highly  compressed  detona¬ 
tion  product  gases.  However,  it  is  rapidly  converted  into  kinetic  energies  of 
gas  and  metal  as  the  case  expands  until,  at  the  moment  of  bursting,  most  of  it 
is  in  this  form.  Since  Gurney  knew  from  experiments  that  radial  fragment 
motion  completely  dominates  axial  components,  he  partitioned  EC  into  radial 

2 

energies  of  metal  and  gas,  integrating  over  1/2pv  for  the  latter.  He  took 
Equation  (15)  above  to  be  true,  using  the  symbol  v  for  vn,  the  case  velocity 

O  D 


7 r 
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at  the  moment  of  bursting  (equal  to  the  gas  velocity  at  the  case),  and  the 
symbol  a  instead  of  rR  for  the  case  radius  at  this  time.  He  then  wrote  his 
Equation  (1): 

2 

EC  =  SsMv  2  >s  vR2  f  B  2nrp  dr  •  (25) 

B  J°  rB 

Here  we  are  using  M  for  the  case  mass  per  unit  length,  a  sum  over  mass  ele¬ 
ments  which  form  a  ring.  Since  his  basic  assumption  requires  EC  to  be 
Independent  of  r0,  the  integral  in  Equation  (25)  must  somehow  be  independent 

of  r  .  Near  the  top  of  page  5  in  his  report,  Gurney  states  that  he  assumed  p 

B  2  2  2 
to  be  constant.  Clearly,  he  took  p  -  C/{ wr  )  =  p  r  /r  ,  a  constant  which 

BooB 

can  be  taken  outside  the  Integral  sign,  putting  rD  in  the  denominator. 

D 

Gurney's  reason  for  assuming  a  linear  dependence  of  v  on  r  now  becomes  clear. 

Equation  (15)  puts  r2/rR2  in  the  integral.  Since  the  integral  of  r^  is  r^/4, 
l\  k 

r  appears  in  the  numerator  to  cancel  r  in  the  denominator  and  make  EC 
a  B  2 

independent  of  the  projectile  size.  The  integral  becomes  CvD  /4  and  a 

D 

solution  for  Vg  gives  us 

vQ  =y]  2E/  [(M/C)  +  1/2],'  (26) 

which  is  Gurney’s  widely  used  cylinder  formula. 

Gurney  was  also  interested  in  small  projectile  warheads  and  grenades  which 
more  closely  resemble  spheres  than  long  cylinders.  His  Equation  (2)  for  a 
cased  spherical  charge  is  analogous  to  his  Equation  (1): 

r  2 

EC  -  S  Mvr2  *  %  Vg2  f  B  4ur  p  j  dr  •  (27) 

B  rB 

o 

where  K  is  now  the  case  mass  and  C  =  p  (4/3^r  J)  is  the  charge  mass  with  r  the 

oo 

spherical  radial  coordinate.  Again  he  took  P  =  3C/(1Jwr  J)  to  be  constant 

5  B  n  5 

which  put  rD  in  the  denominator.  Of  course  the  integral  of  r  is  r  /5  which 

5  B 

puts  r  in  the  numerator,  making  EC  independent  of  projectile  size  in  this 
1:1  2 

geometry  also.  The  integral  becomes  3Cv  /TO,  and  a  solution  for  v  gives  us 

D  B 

Vg  =yj  2E  /  [  (M/C  j  +  3/5] 7  (28) 

for  a  sphere.  Gurney  noted  that  Equations  (26)  and  (28)  make  v  vary  as 

•5  .25  B 

(C/M)  for  small  C/M  and  as  (C/M)*  for  large  C/M  in  agreement  with  the 

obaerva'lons  he  set  out  to  explain  as  simply  as  he  could.  Of  course  both 
equations  tend  to  asymptotic  values  as  (C/M)-»®  (bare  charge). 
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In  his  discussion  in  Section  IV  of  his  report,  Gurney  suggests  a  slightly 
lower  value  than  1/2  in  Equation  (26).  He  also  comments  that  his  equations 
are  intended  to  express  the  fact  that  a  fraction  of  the  chemical  energy 
released  is  converted  into  kinetic  energy,  other  details  being  unimportant  for 
his  purpose.  (His  report  contains  the  word  "important"  instead  of  unimpor¬ 
tant.  This  is  an  obvious  typographical  error  since  the  context  demands  the 
latter.)  In  addition,  he  remarks  that  his  taking  rg  as  the  upper  limit  of 


his  integral  suggests  a  dependence  on  metal  strength.  He  goes  on  to  assert 
that  he  did  this  to  simplify  the  calculation,  since  the  evidence  available  to 
him  at  the  time  gave  no  hint  of  such  a  dependence.  Later  experiments  have 
confirmed  the  fact  that  there  is  a  small  dependence  of  launch  speed  on  metal 


i|2 

strength.  For  example,  Famigliettl  found  that  for  identical  geometries  and 
explosive  fills, the  launch  speed  of  Hadfleld  manganese  steel  fragments  was 
about  15%  greater  than  that  of  stainless  steel  fragments  in  spite  of  the  fact 


that  C/M  was  virtually  the  same  in  both  cases.  This  implies  that  Gurney's  ^2E 
is  not  strictly  a  characteristic  of  the  explosive  alone,  although  it  has 

11-14 

usually  been  treated  as  if  it  were. 


In  Section  V  of  his  report,  Gurney  notes  that  the  value  yf2p=  2.44  x  105 
cm/s  gave  a  good  fit  for  TNT-filled  shells,  over  the  range  0.18  <  M/C  <  16.67, 

17 

which  is  a  bit  larger  than  the  range  examined  by  Taylor  and  Jones.  He  also 
notes  that  this  fitted  value  gives  an  E  value  which  is  only  80$  of  the  energy 
per  unit  mass  which  ought  to  be  released  in  a  TNT  detonation.  This  is,  of 
course,  expected  since  E  is  only  the  contribution  to  the  total  kinetic  energy, 
the  rest  of  the  energy  still  remaining  as  internal  energies  of  gas  and  metal 
after  fragmentation  (neglecting  minor  factors, such  as  light  emission  which 
occurs  later  anyway). 


In  the  last  paragraph  of  his  report,  Gurney  mentions  Taylor’s  paper  and 
correctly  states  that  Taylor  neglected  radial  gas  motion  compared  to  longitu¬ 
dinal.  Most  likely  he  had  not  read  Taylor’s  paper  at  this  time  and  had  this 
on  hearsay  since  he  continues  by  saying  that  Taylor's  results  "were  not 
intended  to  apply  to  a  projectile  from  which  the  end  sprays  are  feeble 
compared  with  the  side  spray."  Gurney  was  obviously  not  aware  that  Taylor 
neglected  v  compared  to  u  in  moving  coordinates  but  found  Vp  >>  Vx  in  labora¬ 
tory  coordinates  (small  angle  <j>  in  Equation  (11)).  Gurney  was  a  physicist 
with  remarkable  insight  who  could  grasp  the  essentials  of  a  problem  without 
much  mathematical  analysis.  For  an  appreciation  of  Gurney,  especially  his 

43 

insights  into  quantum  mechanical  tunneling,  see  Condon.  As  we  have  seen  in 
our  discussion  of  Equations  (14)  and  (15)  above,  Gurney  was  making  an 
excellent  approximation  in  Equation  (15).  In  effect,  he  was  letting  u  =  D  in 
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E.  U.  Condon,  "Tunneling  -  How  It  All  Started,"  American  Journal  of  Physics, 

46,  1978,  p.  319. 
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o 


Equation  (13)  when  he  let  P  =  P_  =  C/(irr2),  which  is  very  nearly  true  for 

bursts  which  occur  for  r/rQ  >  1.1 6  (see  Table  1).  As  Taylor  pointed  out,  the 

gas  density  is  nearly  a  constant  over  an  entire  plane  at  any  distance  x  from  a 
moving  observer,  that  is,  for  any  given  r/rQ.  For  a  sphere  in  which  only 

radial  motion  occurs,  the  analog  of  Equation  (14)  is 


1 

2 


_d 

dr 


(r2pv)  =  0 


(29) 


which  can  be  considered  together  with  mass  conservation 


r  4  3  4  3 

C  =  po  3  ^o  =  p  3  TO  (30) 

to  arrive  at  Equation  (75)  in  this  geometry  also.  If  y  is  the  cartesian  coor¬ 
dinate  normal  to  an  initially  planar  explosive/raetal  interface  and  the  detona¬ 
tion  propagates  along  the  x-axis,  the  analog  of  Equation  (14)  is 

3  3 

"Sy  (pv)  +  (pu)  =  o  ,  (31) 

where  each  term  in  Equation  (31)  vanishes.  Mass  conservation  requires p  uy  s 

p  Dy  where  y  is  the  initial  explosive  sheet  thickness,  so  u  is  not  a 
0  0  o 

function  of  x.  For  u  «  D,  pwl/y  so  v  ~  y  to  make  pv  independent  of  y. 

44 

Thomas,  a  co-worker  of  Gurney,  compared  Gurney's  and  Taylor's  models. 

45 

The  year  before,  he  originated  the  method  of  mentally  dividing  a  shell  or 
bomb  into  sections  by  planes  perpendicular  to  its  longitudinal  axis  in  order 
to  apply  Gurney's  method  of  calculating  launch  speeds,  taking  into  account  the 
curvature  and  variable  wall  thickness  of  real  shells.  His  method  is  still 
used  today  and  has  been  re-invented  several  times,  most  recently  by  Burman  and 

46  44 

Bedford.  In  his  synthesis,  Thomas  began  with  a  brief  description  of 
Gurney's  model,  generalizing  Equations  (26)  and  (26)  to  include  the  planar 
case 

V  =  v0  =V  2E  /  [(M/C)  +  n/(n  +  2)]/  (32) 

where  n  =  1,  2  and  3  for  planar,  cylindrical  and  spherical  symmetries 
respectively.  Following  this,  he  gave  a  brief  account  of  Taylor's  model,  but 
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used  the  exact  curvature  d<£/ds  =  sin<£(d<£/dr)  as  Allison  and  Shriempf  did 
sixteen  years  later.  Thomas  went  on  to  solve  the  Lagrangian  equations  of 
motion  by  a  series  of  successive  approximations,  using  a  parameter  inverse  to 
the  detonation  velocity.  He  obtained  Gurney's  result  as  the  asymptotic  limit 
of  the  zeroth  order  solution  in  which  detonation  is  imagined  to  occur  simul¬ 
taneously  everywhere  in  the  explosive  (D  =  “).  Taylor's  model  emerged  in 
higher  order  solutions.  He  concluded  that  Taylor's  model  should  be  an  excel¬ 
lent  approximation.  Thomas  also  considered  shocks  and  concluded  that  their 
effect  on  the  energy  is  slight.  After  two  or  three  brief  reverberations,  the 
motion  settles  down  to  an  asymptotic  form.  For  very  thin  shells  stepwise 
shock  acceleration  at  very  early  times  has  been  observed  by  Allison  and 

Shriempf, ^  Eden  and  Wright^  have  observed  a  similar  effect  for  thin  plates. 
Theoretical  studies  of  shock  acceleration  have  also  been  carried  out  in  recent 

48  49 

years,  using  the  method  of  characteristics.  *  Neither  of  these  studies 
changes  Thomas'  conclusion  that  there  is  no  need  to  consider  shock  effects 
when  considering  the  motion  of  most  practical  devices. 

Thomas  went  on  to  consider  planar,  cylindrical  and  spherical  cased  charges 
initiated  simultaneously  at  all  points  on  a  central  plane,  on  an  axis  of 
symmetry  or  at  the  center  respectively.  For  the  planar  case,  he  pointed  out 
that  the  problem  is  equivalent  to  Lagrange's  problem  in  Internal  ballistics 

50 

which  was  treated  in  detail  by  Love  and  Pidduck »  This  theme  was  later  taken 

up  by  Jacobs, ^  Sterne^2 later  extended  Gurney's  work  to  explosive  sand¬ 
wiches  as  well  as  cored  cylinders  and  spheres  with  a  fuze  cavity.  He  also 

38 

took  up  Thomas '  model  of  symmetrical  initiation  as  did  Gurney . 
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Let  us  compare  Gurney’s  model  with  Taylor’s.  Gurney's  model  is  more  ele¬ 
mentary  since  It  relies  on  an  energy  partition  at  a  particular  moment,  while 
Taylor  solves  an  equation  of  motion.  With  Gurney’s  model,  we  can  only  estimate 
a  speed  at  a  given  time.  With  Taylor's  model,  we  can  calculate  the  position 
(case  shape)  and  velocity  (speed  and  direction)  of  each  case  element  at  any 
time.  For  many  purposes,  a  Gurney  speed  with  or  without  the  addition  of  a 
Taylor  angle  is  sufficient.  If  metal  trajectories  are  required,  Taylor's 
model  might  be  preferred.  Thomas'  solution  of  the  Lagrangian  equations  could 
be  used,  of  course,  but  this  is  usually  too  complicated  an  approach  for  design 
work  where  the  aim  is  Insight  with  minimum  effort  rather  than  great  precision. 
Sometimes  it  is  useful  or  even  necessary  to  Include  shock  effects,  as  when  one 

54 

metal  plate  is  used  to  accelerate  and  fragment  another. 

More  detailed  modeling  might  Include  microscopic  defects,  crack  propaga- 

55-58 

tion  or  even  molecular  interactions . 


V.  TAYLOR'S  MODEL  SIMPLIFIED  AND  EXTENDED 
A.  Cylinder  Struck  Slde-on 

Since  the  publication  of  Taylor's  papers  by  Batchelor  in  1963,^  his  model 

2 

has  been  more  readily  accessible  to  a  wider  audience  than  Gurney's  report, 
which  was, practically  speaking, an  internal  laboratory  memorandum  which  rela¬ 
tively  few  people  have  actually  read.  In  spite  of  this,  Gurney's  model  has 
been  widely  used  and  extended,  while  Taylor's  model  has  been  virtually  ig¬ 
nored,  except  for  one  equation.  Clearly,  this  is  because  of  the  relative  sim¬ 
plicity  of  Gurney's  model.  In  this  section  we  will  simplify  Taylor's  model  in 
order  to  facilitate  its  use  in  problems  for  which  it  is  appropriate. 
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As  we  mentioned  after  Equation  (25)  above,  Gurney  took  P  to  be  uniform 
inside  the  case  for  any  l  =  r/rQ  value.  In  effect,  he  took 


P  = 


(33) 


where  r  =  r._  (late  in  the  expansion).  Equation  (13)  can  be  written 

O 

P .  CP0D)/t2u)  -  p0/cA)  ,  (jl|) 

where  f  =  u/D  «  1.  Equations  (34)  and  (33)  are  the  same  for  f  =  1.  If  we  put 
Equation  (34)  for  p  in  Equation  (1)  with  y=  3  and  Equation  (34)  for  (Pu) 

In  the  first  term  of  Equation  (20),  we  find 


pu  -  [{BP/)/(Df3)]/£4  =  U3/*4  (35) 

which  defines  U  .  From  Figure  1  for  a  TNT  loading  density  near  that  discussed 

B  10  233  232 

by  Taylor,  we  see  that  B  =  2  x  10  (dyne/cm  )/g/cm  )  =  2(mm4ts)  /(g/cm  )  . 

3  3 

Consequently,  for  P  =  1.51  g/cm  and  D  =  6.38  mm//*s,  UD  =  .7/f  =  1  mm/fis  if 

O  D 

f  =.  0.89*  This  value  of  f  can  be  taken  as  an  effective  value  during  the  ac¬ 
celeration  period  and  enables  us  to  deal  with  UR  as  a  constant  dependent  only 
on  explosive  properties. 

Next  let  us  approximate  u  by  the  formula 

u  =  It  -  Ub(1/£2  +  1A4)  ,  (36) 

where  Um  =  6.914  mm/^s  from  Table  1.  If  UR  =  1  mm/  s,  Equation  (36)  can  be 
used  in  .Equation  (34)  to  find' 

p=  9.634/[6.9l4Ji2  -  (1  +  l/*2)]  ,  (34a) 

3  3 

where  PqD  =  (1.51  g/cm  )  x  (6-38  mm/^s)  =  9.634  (g/cm  )(mm/fis). 

In  Figure  3  we  plot  p  versus  l  from  Table  1  as  the  solid  line  and  Equation 
(34a)  as  the  long  dash  line  which  deviates  from  the  solid  line  by  at  most  15? 
(near  2.  =  i .  1 )  and  Is  indistinguishable  from  it  on  the  scale  shown  for  l  > 

1.5.  A  somewhat  smaller  UR  would  improve  the  agreement  except  near  £  =  1 . 

2  2 

The  short  dash  line  in  Figure  3  plots  P  =  pQ/l  -  1.51  /A  ,  showing  why 

Gurney's  cylinder  formula  based  on  a  radial  energy  partition  can  be  applied 
during  the  later  stages  of  the  expansion  of  devices  initiated  at  one  end.  A 

plot  of  P  =  P^/l2  =  2/22  would  lie  considerably  higher. 

Figure  4  plots  u  versus  l  from  Table  1  as  a  solid  line  and  Equation  (36) 
with  the  values  mentioned  as  the  dashed  line.  Smaller  UR  would  make  the  dash¬ 
ed  line  lie  higher.  The  per  cent  deviation  is  about  the  same  as  in  Figure  3» 
reaching  a  maximum  of  about  15?  near  2  =  1.1.  Disagreement  with  Taylor's  cal¬ 
culations  during  the  early  stages  of  the  motion  is  not  necessarily  bad.  Be¬ 
cause  he  used  a  model  in  which  the  width  of  the  reaction  zone  is  neglected,  we 
expect  his  results  to  be  in  error  at  early  times.  This  point  will  be 
discussed  further  in  Section  C  below. 
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Figure  3.  Gas 
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Figure  4.  Moving  Coordinate  Axial  Gas  Velocity  in  an  Expanding  Cylinder. 


Now  let  us  put  Equations  (35)  and  (36)  In  Equation  (20)  and  require  V  =  0 
for  l:  1.  We  find 


Uw  s  D  +  Ug  ,  (37) 

so  that  Um  can  be  estimated  for  any  explosive  loading  from  D,  pQ,  B  and  f  — 

knowledge  which  Is  available  for  most  cases  of  Interest.  We  may  assume  f  = 
0.89  or  use  Taylor's  f  =  u/D  a3  a  universal  function.  Equation  (20)  becomes 

V  =  P0f’3/2  V(2B/«)(  1-1/j,2/  (38) 

which  does  not  depend  on  D  explicitly  as  In  Taylor’s  model  (only  through 

and  B).  The  parameter  a  =  M/C  appears  In  the  denominator  under  a  square  root 
sign  as  In  Gurney's  formula,  although  there  Is  no  constant  added  to  since 
Taylor  neglects  v  compared  to  u.  This  will  be  discussed  below  In  Section  VI. 

As  l->°°  and  f  «  1,  V  =  PQ  ^/(2B)  fa  and-y2B  plays  a  similar  role  to  yj  2E  'in 

Gurney's  formula.  Here,  however,  there  Is  an  explicit  dependence  on  pQ,a s 

well  as  that  implicit  In  B  and  a  .  It  would  be  desirable  to  tabulate  values 

- T7 

2B  P  for  cases  of  practical  interest.  However,  to  do  a  proper  job,  we 
o  • 

should  make  a  critical  comparison  of  p  versus  p  curves  calculated  by  various 

methods.  Alternatively,  we  can  tabulate  ^2B P  ^ *  from  experimental 
Information. 


2  3  2 

If  we  use  Taylor's  values  of  f  in  Equation  (38)  and  B  =  2(mm/;is)  /(g/cm  ) 
with  a  s  2.67  for  the  standard  cylinder,  we  find  the  values  given  in  column  3 
of  Table  2.  These  are  about  15f  higher  than  the  experimental  values  which 
have  been  repeated  In  colqmn  2.  If  we  let  B  =  1.5  Instead  of  2,  we  find  the 
values  In  column  4  which  exhibit  a  root  mean  square  error  of  0 . 04  mn/ps  rela¬ 
tive  to  the  six  experimental  values.  Here  we  are  lowering  B  by  25%  much  as 
Gurney  lowered  E  from  the  value  he  expected. 

Since  V  »  Vp  =  dr/dt  =  rQdil/dt,  we  can  integrate  Equation  (38)  with  f 
constant  (Ug  constant)  to  find 

Dt»  r  Vf  C«D)/(2Un)  H )*  as  x  (39) 

O  D 

since  v^  =  dx/dt  =  D  cos  <f>  »  D  as  was  mentioned  above  In  connection  with  Equa¬ 
tion  (10).  If  we  put  Equation  (39)  Into  Equation  (4),  it  Is  clear  that  a  case 
element  which  is  seen  by  a  laboratory  observer  to  have  very  little  axial  mo¬ 
tion  (X  w  Xq)  is  said  by  the  moving  observer  to  recede  axially  at  about  the 

detonation  speed.  Of  course,  both  observers  see  the  same  radial  motion.  We 
can  rewrite  Equation  (39)  a3 

r2/r  2  -  x2/[r  2«D/(2Un) J  =  1  ,  (40) 

O  0  D 

which  exhibits  the  approximately  hyperbolic  form  of  the  case  shape  in  moving 
coordinates. 
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Since  V^/V^  =  dX/dr  =  tan  (</>/2)  from  Equation  (11)  and  tan  4>  ~  2  sin  ( <t>/2 ) 

ss  2  tan  (<#>/2),  we  can  integrate  Equation  (38)  with  U0  constant,  using  V  =  2D 

sin  (tf>/2)  from  Equation  (20)  to  find  the  X,  r  path  of  a  case  element  in 
laboratory  coordinates,  namely, 

*  -  x0  =  r„  ta"  W2)"  '  -  1  -  COS-1  Cl/«]  (111) 

which  also  gives  us  X(t)  since  J,(t)  is  known  from  Equation  (39)  •  Accelera¬ 
tions  as  functions  of  time  may  also  be  found.  High  speed  cinematography  or  a 

17 

sequential  series  of  still  photographs  as  used  by  Taylor  snd  Jones  could  be 
used  to  check  some  of  these  relationships. 

Column  5  of  Table  2  gives  a  more  complete  listing  of  Taylor's  values  of 

£  /T7 

ro Ya  than  column  6  of  Table  1,  while  column  6  of  Table  2  give3  values  calcu¬ 
lated  from  Equation  (39)  using  D  =  6.38  ram/^s  and  1)^  =  .75  roa/p-s  corresponding 

to  B  =  1.5  instead  of  2.  Column  7  gives  values  for  the  same  caae  shape 
function  x(  *■)  as  will  be  explained  in  part  C  below. 

Column  8  of  Table  2  gives  experimental  values  of  the  time  t(^)  observed  by 

Lee  and  co-workers^  for  TNT  with  p  =  1.63  g/cm^  («  =  2.47)  in  the  standard 

o 

copper  cylinder.  Columns  9  and  10  give  values  of  t  calculated  from  Equation 
(39)  using  r  =15.3  mm  for  this  test,  D  =  6.38  mm/^s  and  Un  =  1  or  .75  mm/f*s 
o  2  3  2  “ 

corresponding  to  B  =  2  or  1.5  (mm/us)  /(g/cmJ)  .  During  the  later  stages  of 

the  expansion,  the  calculated  values  tend  to  bracket  the  experimental  values 
but  are  higher  in  the  early  stages. 

Column  11  of  Table  2  gives  experimental  values^  for  64/36  RDX/TNT  with 
■a 

P  =  1.717  g/cm  ( «  =  2.34)  in  the  standard  cylinder.  From  Figure  1  we  aee 
0  2  3  2 

that  B  =  2.4  (mm/^s )  /(g/cm  )  for  this  loading.  Column  12  gives  V  calculated 
from  Equation  (38)  for  B  =  2.4,  while  column  13  uses  B  =  .7(2.4)  =  1.68.  This 
illustrates  the  fact  that  Equation  (38)  and  presumably  Equations  (39)— (41)  and 
others  related  to  them  can  be  applied  to  other  explosive  loadings.  Compari¬ 
sons  with  experimental  values  covering  a  larger  range  of  a  values  would  also 
be  desirable,  but  little  time -dependent  data  is  available.  In  Section  VI  we 
will  make  some  comparisons  of  Equation  (38)  with  launch  speeds  for  a  range  of 
a  values. 
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Table  2.  Comparisons  of  Simplified  Taylor  Model  with  Experiment 
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B.  Plate  Struck  Side-On 


Now  let  us  apply  Taylor’s  model  to  a  planar  rather  than  a  cylindrical  ge¬ 
ometry.  The  situation  is  similar  to  that  in  Figure  2  except  that  we  will  use 
the  Cartesian  coordinate  y  instead  of  r,  so  tan  </>  =  dy/dx.  We  will  allow  the 
metal  plates  to  have  different  thicknesses  h1  and  hQ.  We  will  let  the 

thickness  of  the  explosive  sheet  be  W  =  2y_.  The  mass  per  unit  area  of  each 


metal  plate  will  be  designated  by  = 


PMihl  Wlth  1  s  0>1  where  allow  the 


plates  to  have  different  densities  as  well  as  different  thicknesses. 
Equation  (6)  is  now  replaced  by 


M.  D2  -w~ 

1  4 

. 


dx) 

TT3/2 


a  p . 


We  multiply  each  side  of  this  equation  by  2dy/(M1D2)  =  2  tan  if>  dx/CM^D2)  and 
obtain 

2 

2  sin4>  dtj>  =  - =■  pdy 

M.D2  '  (43) 

l 

where  the  left  side  is  readily  integrated  as  in  Equation  (8).  Equations  (10) 
and  (11)  will  have  V  instead  of  Vr  but  are  otherwise  unchanged.  Equation 
(13)  is  replaced  by  y 

PoDyo  =  Puy  ,  (im) 

and  Equation  (1*1)  by  Equation  (31)  which  we  have  already  discussed.  Equations 
(16)  through  (19)  are  unchanged,  so  the  analog  of  Equation  (20)  is 

2  /*y  2p.°y.  r  /  ,\ 


if 


pdy  = 


■  [w  ’  u  '  D] 


2  2 
-  4  sin  (4>/2)  tan  <{> 

=  (V/D)2 


where 


cx.  =  M./(2yoPo)  =  M./C  (46) 

and  C  is  the  explosive  mass  per  unit  area.  We  have  not  repeated  several  steps 
here  since  they  are  the  same  as  before.  The  analog  of  Equation  (23)  is 


-fy  dy/tan 
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Now  let  us  simplify  this  model.  Since  u  is  a  function  of  y  by  Equation 
(44)  as  is  p  =  PQ/(f£)  where  £  =  y/yQ  and  f  s  u/D,  we  can  write  the  analog  of 


Equation  (35) 


where  U_  has  the  same  form  ag  in  Equation  (35). 
B 


(48) 


Now  let  U3  approximate  u  by  the  formula 


u  = 


(49) 


which  is  the  analog  of  Equation  (36).  Here  the  inverse  powers  of  f,  are  n  and 
2n  with  n  =  1,  Just  as  they  were  in  Equation  (36)  with  n  =  2.  The  spherical 
analog  might  consist  of  a  disk-shaped  detonation  front  rotating  about  a 
radius,  s  configuration  very  difficult  to  realize  experimentally.  If  we  use 
Equation  (49)  in  Equation  (44),  we  find 


p  •  p0°/(y  -  u,  (i  *  I)]  •  (so) 

If  we  use  p  D  s  9.634(g/cm^)(mm/p.s) ,  Ura  =  6.914  mm/p.s  and  U  =  1  mm/ps,  we 

obtain  the  decreasing  long  dash  line  in  Figure  5  which  diminishes  more  slowly 

than  the  analogous  curve  in  Figure  3*  The  short  dash  line  just  above  it  cor¬ 
responds  to  p  =  p^/i  s  2/£  for  the  head-on  case,  while  the  dotted  line  below 


it  is  for  p  =  pQ/£  =  1.5 1/A ,  a  Gurney-type  assumption.  The  Increasing  da3hed 

line  plots  Equation  (49)  for  the  same  Uo,  and  U0.  This  line  increases  more 
slowly  than  its  analog  in  Figure  4. 


Next,  let  us  put  Equations  (48)  and  (49)  in  Equation  (45)  and  require  V  = 
0  for  £  s  1,  giving  Equation  (37)  again.  Equation  (45)  becomes 


v  ‘  v3/V(.f)  l1 '  i)  <5° 

which  is  the  analog  of  Equation  (38)  with  n  =  1  replacing  n  =  2  as  the  coeffi¬ 
cient  of  B  and  the  inverse  power  of  £.  As  £-*»  and  f  ^  1  we  can  compare  Equa¬ 
tion  (51)  with  Equation  (32)  above  with  n  =  1.  For  two  equal  mass  metal 
plates  in  a  symmetric  sandwich,  a  =  2  a1  =  2M/C  in  Equation  (46),  so  Equation 

(51)  becomes  V  =  P  >/2BA*  in  this  limit.  We  will  compare  these  formulas  for 

a  range  of  ct  values  in  Section  VII  below. 


Since  V  aVy  ;  dy/dt  =  yQdt/dt,  we  can  Integrate  Equation  (51)  with  f 
constant  to  find  the  analog  of  Equation  (39) 


Dt 


w  yQ-y  -^j-  |i>/i  -l/l  +  h  In  [lyjl  -  l/l  +  J,-Ji]-%*nJsj  (52) 


«  x 


» 
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u  (mm/^is) 
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Figure  5.  Gas  Density  and  x  Velocity  Component  for  an  Expanding  Sheet. 


Clearly  the  case  shape  (y  versus  x)  is  not  approximately  hyperbolic  as  it  was 
for  the  cylinder.  Since  V  /V  =  dX/dy  =  tan(<f>/2)  by  the  analog  of  Equation 

x<v  y 

(11)  and  tan  <i>  ^2  sin(<l>/2)  ^  2  tan(4'/2),  we  can  Integrate  Equation  (45)  with 
(or  f)  constant  to  find  the  X,  y  path  of  a  plate  element  in  laboratory 

coordinates, 

x  -  xo  =  y0VV(4aiD^  {W1  '  1/1  -  h  inW  -  +  *  -  h]  +  %  in  ( 

which  is  the  analog  of  Equation  (41).  Use  of  Equation  (52)  gives  X(t). 
Accelerations  may  also  be  found. 

C.  Expansion  in  the  Reaction  Zone 


As  we  noted  above,  Taylor  used  a  model  which  assumes  a  zero  width  for  the 

21 

reaction  zone.  H.  Jones  has  pointed  out  that  case  expansion  occurs  inside 
the  reaction  zone,  leading  to  a  slight  decrease  in  detonation  velocity  com¬ 
pared  to  the  idealized  Chapman-Jouget  value.  He  also  derived  an  expression 
for  the  shape  of  an  expanding  cylinder  valid  inside  the  reaction  zone: 


VI 


-1, 


—  —  =  2  cosh  (i) 
a  r 

o 


(54) 


for  1  <  i  =  r/r  <  r„/r 
—  o  —  1  o 


Here  we  have  used  a  as  defined  in  Equation  (21) 

above  instead  of  Jones'  mass  per  unit  area  a  =  M/(2nr  ).  At  the  rear  of  the 

o 

reaction  zone  we  can  rewrite  Equation  (54)  as 


h  *  Vro  ■  C05h  S;VI=  1  *  *  ■  (55) 

Jones  gave  a  numerical  example  in  which  he  took  x^/rQ  =  2/a  =  1/2,  assuming 

that  the  width  of  the  reaction  zone,  x^,  might  be  as  large  as  half  the  charge 

radius  for  a  =  4.  Thus,X  =  0.015  in  Equation  (55)  so  a 1  is  only  slightly 

larger  than  unity.  This  would  lead  to  a  6.25$  reduction  in  detonation  veloc¬ 
ity  compared  to  the  idealized  value.  For  a  bare  charge  (the  extreme  departure 
from  a  rigid  tube), Jones  estimates  a  9.6$  reduction  for  the  same  reaction 
zone.  For  solid  or  liquid  explosives, the  idealization  of  a  rigid  tube  cannot 
be  realized  experimentally  since  at  least  the  inner  radius  of  even  a  very 
heavy-walled  tube  of  metal  will  expand  somewhat  under  the  pressures  produced 
by  condensed  explosives  (unlike  a  shock  tube  filled  with  a  detonable  gas 
mixture) . 


Column  7  in  Table  2  exhibits  Equation  (54)  as  a  function  of  i  and  compares 
it  with  Taylor's  values  in  column  5, as  well  as  the  present  model  values  in 
column  6.  According  to  Jones  there  is  less  expansion  near  the  front  of  the 
reaction  zone  than  Taylor  estimates  (or  a  given  amount  of  expansion  occurs 
farther  behind  the  detonation  front).  However,  the  two  curves  agree  rather 
closely  even  in  the  reaction  zone  (say  1  <  a  <  1.03)  and  deviate  appreciably 
only  well  behind  the  reaction  zone  where  Jones'  formula  is  no  longer 
applicable.  This  is  illustrated  in  Figure  6. 
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VI.  THOMAS'  HEAD-ON  MODEL  SIMPLIFIED 


44 

Thomas  developed  his  own  model  for  plane,  cylindrical  and  spherical 
cased  charges  initiated  on  a  symmetry  plane,  axis  or  point.  In  these  cases  the 
detonation  front  will  strike  all  points  of  the  metal  surface  simultaneously, 
"head-on"  (propagation  vector  parallel  to  surface  normal).  If  dm-  is  the  mass 
of  product  gas  in  a  volume  element  of  thickness  dr  (per  unit  area  for  the 
plane,  per  unit  length  and  angle  for  the  cylinder  and  per  unit  solid  angle  for 
the  sphere),  then  this  specific  mass  element  Is 


,  -  „  n-1, 

dm  =  pr  dr  , 

where  n  =  1,  2  and  3  for  plane,  cylinder  and  sphere  respectively.  The 
equation  of  gas  motion  ia 


32r  _  _  I  iE  =  .  -*>-1 
9t2  P  3l> 


3r  3p  n-1  3p 
3m'  T  3m' 


(56) 


(57) 


Thomas  used  Equation  (1)  for  p  but  kept  y  arbitrary  and  concentrated  on  cylin¬ 
drical  symmetry.  He  was  able  to  obtain  a  first  integral  but  could  not  com¬ 
pletely  solve  the  problem.  He  assumed  an  initial  uniform  distribution  of  mass 
and  outward  velocity  when  the  detonation  wave  reached  the  metal  surface  with 
boundary  conditions  m'  =  0  for  r  =  0  and 


M 


r”"1  P 


(58) 


for  m'  =  C,  the  total  specific  charge  mass,  at  the  moving  metal  surface.  As 
before,  M  is  the  specific  metal  mass.  After  a  brief  period  of  shock  reverber¬ 
ation,  Thomas  argued,  the  solution  ahould  approach  a  form  in  which  the  method 
of  separation  of  variables  can  be  applied.  Letting  r  be  a  product  of  a 
function  of  time  and  a  function  of  mass,  namely. 


r  =»  £.(t)g(m') 


Equation  (57)  separates  into 


and 


(59) 


(60) 


(61) 


2 

where  we  have  used  v  for  the  separation  constant  instead  of  letting  it  be 
unity  as  Thomas  did.  Dimensionally  speaking,  v  is  an  inverse  relaxation  time. 
The  boundary  conditions  become  g(0)  =  0  and 


2  n  n-1 
M\>  g  =  Bg 


(62) 
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for  m"  =  C  where  we  have  used  Equations  (59)  and  (60)  In  Equation  (58). 

If  we  multiply  Equation  (60)  by  d£  =  £dt  and  Integrate,  we  obtain  the 
first  integral  of  Equation  (60) 


£  =  v 


[1  -  )Tn(Y-1)]  -  V/ro  , 


(63) 


where  we  have  taken  £  =  0  for  £  =  1.  The  second  equality  comes  from  taking 
g(C)  =  rQ  in  Equation  (59)  so  V  =  r  =  rQ£  for  the  case  velocity.  If  Y  =  3  and 

£  ->  00  in  Equation  (63),  V  =  r^v/y/n.  As  £  ■*  00  and  f  ^  1,  Equation  (38)  for  a 

cylinder  and  Equation  (51)  for  a  symmetric  sandwich  (a  =  2a  )  both  give  V  = 

PQ  v/(2B) /a  .  If  we  assume  that  the  head-on  and  side-on  formulas  are  equal 

under  these  conditions  (same  launch  velocity), 


(po/ro)  y[(2B)/a 

and  Equation  (63)  becomes  for  Y  =  3 

V  =  P0^(2B/oO(l  -  r2n) 


(64) 


(65) 


Requiring  approximate  equality  for  large  £  corresponds  to  the  approximate 
equalities  of  the  gas  densities  Illustrated  in  Figure  3  as  £  Increases.  For  a 
cylinder  (n  =  2)  the  head-on  velocity  in  Equation  (65)  may  be  expressed  as 


V  =  Po^(2B/oO(l  +  JT2)(1  -  Z  2) 
and  compared  with  Equation  (38)  for  the  side-on  velocity 


(66) 


V  =  P0^l(2B/a)f"3(l  -  JT2)  (38) 

In  both  cases  V  =  0  for  £  =  1.  They  differ  by  the  factors  Vi  +  £~2  and 
-3/2 

f  ,  where  f  is  given  by  u  in  Taylor's  Table  I  or  Equation  (36)  divided  by 
D.  Similarly  for  a  symmetric  sandwich  (n  =  1),the  head-on  velocity  in  Equation 
(65)  is 


V 


(2B/oO  (1  +  5-"1)  (1  -  2."1) 


(67) 


while  the  side-on  velocity  in  Equation  (51)  is 


V  = 


(2B/ct)f  3(1  -  £-1) 


(51) 


These  formulas  differ  by  the  factors  \l  +  £_1  and  f~^2,  where  f  is  given  by 
Equation  (49)  divided  by  D. 
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Thomas  did  not  let  y  =  3  or  make  such  interpretations  of  £  and  g  and  of 
course  set  v  -  1.  However,  he  did  write  down  the  first  integral  of  Equation 
(61)  for  arbitrary  n  and  y,  namely,  with  K  constant, 

,  2  2  yB  I  n-1  dg  \  "(Y 

h  v  8  +  (T^iy  arj  =  K-  (6, 

For  the  special  case  of  a  cylinder  (n  =  2),  he  also  found  a  second  integral 


m'  =  Y  -  1  [k  Y  '  1  -  (K  -  ^  v2g2)  Y  ’  1  ]  .  (69) 

However,  for  n  =  2  and  arbitrary  y,  he  could  not  integrate  Equation  (63) 
exactly.  Instead,  he  integrated  numerically  with  y  near  3  as  a  parameter. 
Finally  he  found  the  kinetic  energy  of  the  gas  in  the  expanding  cylinder  and 
added  this  to  the  kinetic  energy  of  the  case,  equating  the  sum  to  EC  in  order 
to  compare  with  Gurney's  model.  In  the  limit  £  *>  =>  he  wrote  a  formula  for  V 

which  depended  on  y,  as  well  as  \^2E  and  a.  Even  in  this  limit  his  formula  for 
arbitrary  y  is  rather  complicated.  However,  for  y  =  3  it  simplifies  to 

V  =  V  U2E)/(2a/5)]  [(1  +  l/«)2/3  -  1  ]  /  [(1  +  l/a)5/3  .  i]  (70) 


which  gives  almost  the  same  launch  speed  as  Gurney's  much  simpler  cylinder 
formula, 


V  =  V(2E)/  O  +  h]  .  (26) 

These  two  formulas  are  compared  as  curves  (a)  and  (b)  in  Figure  7.  As  a  ■*  0, 

Thomas'  V  while  Gurney's  V  -*■  V4E.  Both  formulas  approach  zero  as  a  ■*  <*> 

and  are  practically  indistinguishable  for  a  >  0.1.  Of  course  V  will  reach 
zero  for  large  finite  a  when  the  case  fails  to  break  and  the  model  has  ceased 
to  apply. 

Equation  (38)  for  l  -*  f  ^  1  gives  V  =  oq  y/2B/a  which  approaches  zero  as 

a  -*  ®  but  Increases  without  limit  as  n  ■*  0  since  Taylor's  model  takes  no  ac¬ 
count  of  the  gas  energy  and  assumes  a  "sufficiently  heavy"  case  (a  >  0).  One 
way  to  compare  this  form  of  V  with  Gurney's  formula  might  be  to  assume  that 
some  fraction,  say  0.8»of  the  available  energy  appears  as  case  kinetic  energy 
with  the  rest  appearing  mainly  as  gas  energy  as  ^  Then 

%  M M2  =  h  M(2BpQ2/ct)  =  0.8  (EC),  (71) 

and  the  limit  form  of  Equation  (38)  becomes 

V  =  \/(2E)/(ff/0.8)  ,  (72) 

which  is  plotted  as  curve  (c)  in  Figure  7.  If  we  wish  to  avoid  the  large 
divergence  near  a  =  0,  we  might  add  a  constant  like  n/(n  +  2)  from  Equation 
(32)  to  a  in  Equations  (64)  -  (67),  (38)  and  (51).  Then  for  n  =  2, the  limit 
form  of  Equation  (38)  becomes 
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Figure  7.  Comparison  of  Launch  Velocity  Formulas  Versus a  . 
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If  we  require  Equation  (73)  to  reduoe  to  Gurney’s  formula  asnr-»D,  then  V  = 


n/4"e  or 


BPq2  =  E  ,  (74) 

completing  the  analogy  between  B  and  E.  If  we  add  1/3  to  a  for  n  =  1  and 
consider  the  limit  form  (£•+•  °»  f  %  1),  Equation  (51)  becomes  for  =  a/2 

V  =  ^2Bpo2/(a  +1/3)  =  ^2E/(a  +  1/3)  .  (75) 

by  Equation  (74).  This  is  the  same  as  the  Gurney -Thomas  formula  for  a 
symmetric  sandwich.  Equation  (32),  and  is  also  plotted  in  Figure  7.  Equation 
(32)  with  n  s  3  for  a  sphere  might  also  be  plotted  but  has  been  omitted  to 
avoid  further  confusing  the  figure.  For  typical  values  of  a  which  generally 
exceed  unity,  all  of  these  formulas  predict  similar  launch  velocities 
independent  of  geometry  (plane,  cylinder  or  sphere)  or  direction  of  travel  of 
the  detonation  wave  (head-on  or  side-on).  They  differ  mainly  for  small 
None  of  them  are  really  intended  for  a  bare  charge  (a  =  0)  or  for  very  large 
values  of  a  =  m/C.  A  finite  value  for  V  in  the  limit  a  -*■  0  might  be  thought 
of  as  corresponding  to  the  case  in  which  a  finite  mass  of  metal  Is  acceler¬ 
ated  by  an  unlimited  amount  of  explosive.  As  is  well  known,  beyond  a  certain 
point  the  addition  of  more  explosive  does  almost  nothing  to  Increase  launch 
velocity,  even  in  a  vacuum  where  air  resistance  is  not  a  limiting  factor. 
Similarly,  if  the  limit  a  -+■  00  is  approaohed  by  decreasing  the  amount  of 
explosive,  we  reach  a  situation  in  which  a  detonation  wave  will  not  propagate 
because  the  dimensions  are  too  small.  Again  the  model  does  not  apply. 


Thomas  did  not  consider  the  case  y  *  3*  nor  did  he  pursue  a  solution  for 
planar  symmetry  (n  =  1).  However,  this  Is  a  case  for  which  his  formulation 
of  the  problem  may  be  solved  exactly.  For  y  =  3»  n  =  1,  the  Integral  of 
Equation  (63)  is 


£ 


Vi  *  »v . 


(76) 


while  the  Integral  of  Equation  (68)  is 


(77) 


Since  Equation  (77)  gives  g(m')  implicitly  while  Equation  (76)  gives  £(t)  ex¬ 
plicitly,  the  problem  of  finding  r  in  Equation  (59)  is  solved.  The  constant 
K  can  be  found  from  Equation  (68)  by  using  Y  =  3,  n  =  1  and  the  initial 
values  g  =  rQ  and  dm'/dg  =  p^»  the  density  at  the  rear  of  the  reaction  zone. 

We  may  use  Equation  (76)  and  r  =  r  2.  for  £  >  1  in  Equation  (58)  with  n  =  1 
to  find 

MrQv2/£3  =  p  =  Bp3  =  B(Pl/£)3  .  (78) 
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Equation  (64)  with  n  =  1  and  <>  =  2M/C  =  2M/(2roPo)  for  a  symmetric  sandwich 
may  then  be  used  in  Equation  (78)  to  find 

(Pj/  P0>  -  21/3  ■  1-26  (79) 

close  to  the  value  4/3  =  1.33  suggested  by  H.  Jones. 

VII.  STERNE'S  EXTENSIONS  OF  EARLIER  WORK 

52 

In  1947,  Sterne  considered  methods  of  extending  the  ideas  of  Gurney  and 
Thomas  to  asymmetric  sandwiches,  Including  the  extreme  case  of  a  single  plate 
(sometimes  called  an  open-face  sandwich).  He  also  considered  cylinders  with 
solid  metal  cores,  representing  fuzing  devices.  For  two  plates  with  specific 
masses  (per  unit  area)  and  MQ  separated  by  a  sheet  of  explosive  of  specific 

mass  C,  he  wrote  for  the  radial  energy  partition  at  launch 

EC  =  1/2(M1V12  +  MqVo2)  +  (l/6)C(Vt2  +  V02  +  V1V0)  ;  (80) 

and  for  the  conservation  of  momentum  (initially  zero), 

0  =  M,V  +  M  V  +  (172)0(7,  +  V  )  ,  (81) 

110  0  10 

where  >  0  and  Vq  <  0  in  the  chosen  coordinate  system.  Sterne  mentioned 

that  he  was  following  Gurney’s  procedure,  assuming  a  linear  variation  of  gas 
velocity  v  with  distance  and  uniformity  of  P  at  all  points  between  the  plates 
at  launch  time.  He  did  not  bother  to  derive  Equations  (80)  and  (81),  probably 
thinking  it  to  be  sufficiently  obvious.  We  will  give  a  derivation  below  as 
part  of  a  more  general  treatment. 

If  we  eliminate  VQ  between  Equations  (80)  and  (81)  we  obtain  the  asymmet¬ 
ric  sandwich  formula  of  Sterne: 

V,  =  J(2E)/t<*  +  a  q2  +  1/3  (1  -  q  +  q2)]  (82) 

It  1  o 

with  V  =  -qV.  and  q  =  (1  +2  «, )/( 1  +  2a  )  with  «  =  M./C  and  «  =  M  /C.  If 

01  l  olloo 

M1  =  MQ,  =  «Q  =  a/2  and  q  s  1,  then 


V1  =  ^f(2E)/(  «+1/3)  =  -VQ  (83) 

for  the  symmetric  sandwich,  the  same  as  Equation  (32)  with  n  =  1.  If  a  =0, 

_ _ _  .  o 

V1  =  V(2E)/[  1/3  (1  +  5^  +  4  «1 2 )  ]  (84) 

for  a  single  plate.  In  this  limit  the  model  predicts  that  V  =  -  (1  +  2«^)V^, 

so  IV  I  >  V.  as  M  ->  o.  We  can  interpret  V  in  the  limit  H  -+  0  as  the 
o  1  o  K  o  o 

velocity  of  the  free  surface  of  the  gas. 
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Twenty  years  later,  Henry  cited  Sterne’s  report  and  repeated  his 
results.  Unlike  Gurney’s  report,  Sterne's  was  widely  circulated  and  had 
Hughes  Aircraft  among  its  recipients.  Henry's  report  has  been  frequently 
cited  in  the  literature  while  Sterne's  has  been  forgotten.  Probably  Henry's 
report  has  been  the  main  instrument  through  which  most  interested  parties  have 
become  acquainted  with  Gurney's  and  Thomas'  work,  since  almost  no  one  cites 
Sterne  or  Thomas.  It  also  seems  likely  that  Henry  never  read  Gurney's  report 
since  he  does  not  seem  to  have  a  clear  understanding  of  Gurney's  method.  In 
his  own  derivation  of  Equation  (84),  Henry  arrives  at  the  more  cumbersome  but 
equivalent  form 

V1  =  y<2E )/U  +  1/6{1  +  (1  +  2«)3}/d  +  «)]  .  (85) 

The  equivalence  is  easily  seen  by  expanding  the  cubic,  factoring  2  from  the 
curly  brackets,  dividing  by  (1  + «  )  and  collecting  terms.  Similar  considera¬ 
tions  hold  for  Henry's  form  of  Equation  (82)  which  appears  again  more  recently 

9 

in  the  paper  by  G.  Jones  and  co-workers. 
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About  five  years  after  Henry's  report  was  published,  Defourneaux 

2 

repeated  Equation  (84)  as  his  Equation  (18),  citing  Gurney's  report,  but  not 
Sterne's.  Of  course  Gurney's  report  does  not  contain  Equation  (84),  although 
other  references  did  by  this  time.  Defourneaux  pointed  out  that  the  kinetic 

2 

energy  per  unit  explosive  mass  for  a  single  plate,  namely  ^>^  =  (1/2)«^V^  where 

Vj  is  given  by  Equation  (84),  has  a  maximum  (equal  to  E/3)  for  a =  M^/C  = 

1/2.  He  did  not  pursue  this,  but  It  is  interesting  to  investigate  the 
asymmetric  sandwich  as  well.  We  wish  to  optimize  <f>^,  using  Equation  (82)  for 

V^,  subject  to  the  constraints  a ^  >  0,  aQ  >  0.  If  we  let  dty dao  =  0,  we  find 

a,  +  a  +  6a.a  =  0.5  ,  (86) 

1  o  1  o  ' 

which  is  plotted  as  the  dashed  curve  in  Figure  8  and  gives  the  locus  of  the 
inflection  points  in  the  curves  ^(a,j)  for  «o  <  0.5.  Values  of  <f>  in  units  of 

E,  namely,  <f>^/ E,  are  given  in  parentheses  at  some  points  in  the  figure. 

Clearly  <f> ^  increases  with  aQ.  If  we  let  3</>1/3  =  0,  we  find 

otj  =  0.5  ^|tl  +  5a^  +  4ct^)/(l  +  3ct^)  >  (87) 

which  is  plotted  as  the  solid  curve  in  Figure  8  and  gives  a  maximum  <f> .(“.)  for 

o 

any  aQ.  The  straight  line  at  45  is  for  aQ  =  (symmetric  sandwich)  while 
the  horizontal  straight  line  is  for  a 1  =  0.5.  As  a 1  and  aQ  increase  without 


®I.G.  Henry,  "The  Gurney  Formula  and  Related  Approximations  for  the 
High-Explosive  Detonation  Products,"  Hughes  Aircraft  Report  PUB-189*  1967. 

Defourneaux,  "Transferts  d'Energie  Dans  les  Combustions  et  Detonations 
avec Confinement, "  Astronautics  Acta.  17,  1972,  p.  609. 
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Asymmetric  Sandwich. 


limit,  <f>^/E  approaches  unity.  This  incorrect  prediction  says  that  we  achieve 

the  greatest  possible  launch  energy  per  unit  explosive  mass  when  we  sandwich 
the  explosive  between  two  infinitely  heavy  plates.  The  model  fails  for  large 
since  it  does  not  account  for  energy  absorbed  as  heat  by  relatively 

immobile  masses.  This  will  be  discussed  further  below  after  Equation  (104). 
For  the  present,  it  is  sufficient  to  note  that  errors  should  be  less  than  10% 
if  <  7.  For  typical  the  model  3hould  correctly  predict  that  the 

addition  of  a  backing  plate  (aQ  >0)  will  increase  the  specific  launch  energy 

of  plate  1.  Hoskin  and  Lambourn49  arrived  at  a  similar  conclusion  by  using 
the  method  of  reflected  shock  characteristics. 

In  units  of  the  Gurney  velocity, %/5e,  the  launch  velocity  for  a,  =  0.5,  a 

1  o 

=  0  is  V^/v2E  =  2/3  =  0.816,  while  the  optimized  launch  energy  is  </>^/E  = 

1/3  =  0.333.  At  «1  =  =  0.729  where  the  optimized  sandwich  is  symmetric, 

V,/y/2E  =  0.748  (8$  less).  However,  ^/E  =  0.408  (22.5$  greater).  The 
combined  energy  of  both  plates  is  of  course  still  larger.  For  the  optimal 
values  =  0.98,  «o  =  2,  we  find  V^V^T  =  0.72  (down  28$),  but  *  /E  =  0.507 
(up  52$), 
and  so  on. 
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It  is  interesting  to  note  that  Sewell  assumed  an  equivalence  between 
warhead  effectiveness  and  launch  kinetic  energy  (or  momentum)  per  unit 
explosive  mass  and  optimized  the  design  of  cylindrical  and  spherical 
configurations,  using  Gurney's  formulas.  Of  course  these  geometries  are 
necessarily  symmetrical.  His  cylinder  calculations  were  later  refined  by 

Zulkoskl  who  took  account  of  finite  end  effects. 

<5 

G.  Jones  and  co-workers,  who  recently  rederived  Equations  (82)-(84), 

referring  to  Henry's  report60  but  not  to  Sterne's,^2  also  derived  a  formula 

for  the  velocity  of  a  plate  in  a  symmetric  sandwich  based  on  Equation  (83)  and 
Equation  (1): 


V1  = 


1  2E  1*1 

n- (Ywl) (l+l/3a)l 

'  (a  +  1/3)  L1 

-  K  J 

(88) 


If  V  =  3  in  Equation  (88)  and  3 «  >>  1,  this  formula  reduces  to  Equation  (65) 
2 

with  n  =  1,  BpQ  =  E  from  Equation  (74).  For  an  asymmetric  sandwich  they 
obtained 


R.G.S.  Sewell,  "Fixed-weight  and  Fixed-volume  Constraints  on  Optimum 
Charge-to-Metal  Ratios  in  Warhead  Design,"  NAVWEPS  Report  8471,  NOTS  TP 
3430,  China  Lake,  CA,  1964. 

^T.  Zulkoski,  "Development  of  Optimum  Theoretical  Warhead  Design  Criteria," 
Naval  Weapons  Center  TP  5892,  1976. 
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(89) 


V1  ~  +  "  <l]  ~(Y  "  lHq  +  1)Q/alj 

with 

Q  =  ai  +  a0q2  +  1/3(1  -  q  +  q2)  •  (90) 

For  a  single  plate  or  open-face  sandwich  («o  =  0),  q  =  1  +  2 and  Q  =  1/3(1 

+  5 +  4aj2)  in  Equation  (89).  Of  course  y  should  be  3  as  we  have  shown. 

Sterne  derived  similar  formulas  for  V(£)  in  the  case  of  sandwiches.  In  the 
same  report,  he  discussed  "cored”  cylinders,  that  is,  a  Gurney  cylinder  with  a 
solid  metal  core  or  fuze  cavity.  This  will  be  discussed  further  in  the  next 
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sections,  as  will  his  treatment  of  a  "cored"  sphere. 

fill  9 

Chanteret  follows  a  method  similar  to  that  of  Jones  and  co-workers,  but 

treats  the  "side-on"  acceleration  of  plates  and  cylinders  explicitly  in  his 

extension  of  the  Gurney  model.  For  symmetric  cases  he  a^ain  arrives  at 

equation  (32),  the  generalization  first  made  by  Thomas-  However,  his 
expression  for  E  is  the  square  brackets  in  equation  (20)  multipled  by  D.  His 
notation  differs  from  ours  since  his  u  is  our  U,  defined  in  equation  (5) 
above.  His  calculation  of  V  still  involves  finite  time  or  space  steps,  so  it  is 
not  completely  analytical,  but  numerical.  For  an  asymmetric  sandwich  he 
employs  the  usual  technique  of  a  zero  velocity  plane  in  the  explosive.  He 
also  applies  this  technique  to  obtain  a  Gurney  type  formula  for  the  launch 
velocity  of  a  hollow  cylinder,  namely,  a  tube  of  explosive  between  two  tubes 
of  metal,  which  he  looks  upon  as  an  asymmetric  sandwich  rolled  about  an  axis 
above  or  below  the  sandwich  parallel  to  an  edge  until  two  edges  meet.  His 
expression  for  the  location  of  the  zero  velocity  cylindrical  surface  In  the 
explosive  involves  a  cubic  equation  with  the  Chapraan-Jouget  density  of  the 
explosive  products  as  one  of  the  parameters.  His  expression  for  E  still 
Involves  finite  space  or  time  steps,  except,  of  course,  in  the  limit 

In  1965  Hoskln  and  co-workers  presented  experimental  data  for  single 
plate  launch  velocities  in  terms  of  C/M,  the  inverse  of  our  a  =  M/C.  They 
compared  their  data  with  their  version  of  Equation  (84),  namely 

V  =  >/(2E>/(  1/3(1  +  4«x  +  =  \/6E/(1  +  2«t)  ,  (91) 


6lJP.Y.  Chanteret,  "An  Analytical  Model  for  Metal  Acceleration  by  Grazing  . 
Detonation,"  Seventh  International  Symposium  on  Ballistics,  The  Hague-The 
Netherlands,  1993,  p.  51 5. 

65N.E.  Hoskin,  J.W.S.  Allan,  W.A.  Bailey,  J.W.  Lethaby,  and  I.C.  Skidmore, 

"The  Motion  of  Plates  and  Cylinders  Driven  by  Detonation  Waves  at  Tangential 
Incidence,"  Forth  International  Symposium  on  Detonation,  ONR  ACR-126,  1965, 
p.  14, 
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a  simplification  achieved  by  using  4a  instead  of  5 «  in  the  denominator  with 
little  variation  in  «-dependence .  Their  analysis  was  based  on  a  method  of 

6  6 

shock  characteristics  previously  discussed  by  Hill  and  Pack.  Related 
theoretical  work  appears  in  References  48  and  49  mentioned  above.  Equation 
(91)  is  for  the  side-on  case  while  Equation  (84)  is  head-on  since  it  is  based 
on  an  energy  partition  for  motion  normal  to  the  initial  plane  of  the  plate 
surface.  Again  we  see  that  launch  velocity  formulas  (Z  -*■  °°),  in  this  case  for 
single  plates,  do  not  depend  strongly  on  the  direction  of  motion  of  the 
detonation  wave. 

Figure  9  shows  Hoskln’s  data  as  the  solid  line.  Here  we  are  taking  his 
normal  velocity  component  and  to  be  approximately  the  same  as  we  did  with 

similar  cylinder  data.  Our  a 1  is  the  inverse  of  his.  Since  he  gave  no 

numerical  values,  only  a  graph,  the  values  we  show  for  large  <v^  are  somewhat 

uncertain , and  there  may  be  better  agreement  than  shown  between  model  and 
experiment.  Equation  (91)  is  shown  as  the  dashed  line,  using  \/2E  =  2.575  imd/ms 

=  4.46/%/Tas  suggested  by  Hoskins.  On  the  scale  of  the  figure,  it  is 

indistinguishable  from  Equation  (84).  Of  course  Equation  (89)  with  Z  -*■  °°,y  = 

3  and  =  0  is  the  same  as  Equation  (84). 

p 

It  is  interesting  to  note  that  =  (1/2)a,V_  has  a  maximum  for  a  =  1/2 

111  1 

whether  we  use  4  or  5  as  the  coefficient  of  a ^  in  the  denominator  of  Equations 
(91)  or  (84). 


VIII.  SIMPLE  GEOMETRIES  WITH  AN  ARBITRARY  NUMBER  OF  LAYERS 
A.  Introduction 

Henry^0  also  repeated  Sterne’s  modeling  of  a  cored  cylinder"^  and 
mentioned  a  cored  sphere  which  Sterne  also  modeled.  Neither  author  mentioned 
the  planar  analog,  a  sandwich  with  a  metal  center,  which  implies  a  minimum  of 
five  layers  instead  of  three.  Henry  did,  however,  discuss  the  problem  of  a 
cylindrical  configuration  with  more  than  two  or  three  alternating  tubes 
(cylindrical  layers)  of  metal  and  explosive,  an  arrangement  which  he  called  a 
jelly  roll.  If  we  ignore  the  lack  of  layer  closure  in  the  pastry  analog,  this 
is  a  descriptive  name.  The  spherical  analog  might  be  called  an  onion,  if  we 
Ignore  the  fact  that  the  alternate  layers  are  of  different  materials. 

Finally,  the  Dagwood  seems  to  be  a  suitable  name  for  a  multilayered  planar 
sandwich. 

Henryk  wrote  down  a  very  simple  model  for  a  many-layered  jelly  roll  in 
which  he  envisioned  all  metal  layers  to  be  of  the  same  composition  and 


CC 

°R.  Hill  and  D.C.  Pack,  "An  Investigation,  by  the  Method  of  Characteristics, 
of  the  Lateral  Expansion  of  the  Gases  behind  a  Detonating  Slab  of 
Explosive,"  Proceedings  of  the  Royal  Society  of  London,  19],  1947,  p.  524. 
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Figure  9.  Experimental  Single  Plate  Launch  Speeds 
Compared  with  Model  Predictions. 
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thickness  with  metal  on  the  exterior.  He  also  envisioned  all  explosive  layers 

p 

of  the  same  composition  and  thickness.  He  then  wrote  (C  +  M)/(irR  )  for  the 
average  initial  density  of  the  array,  where  C  and  M  have  their  usual  meanings 
for  a  cylinder  and  R  is  the  outer  radius  before  detonation.  His  energy 
balance  became 

E C  -  hf  J"CC  +  M)/(ttR2)1  [v2  )  1  27Trdr 

Jo  L  \R/  J  (92) 

*  1/4  (C  +  M)  V2 

Here  he  used  the  initial  radius  instead  of  the  radius  at  bursting,  apparently 
not  having  Gurney's  report  to  aid  his  understanding  of  the  model.  Equation 
(92)  gives 

V  =  ^(4E)/(1  +  a)  ,  (93) 

where  ®  =  M/C  for  the  entire  device.  He  compared  Equation  (93)  graphically 
with  Gurney's  simple  cylinder  formula,  Equation  (26) 

VG  =  ^(4E)/(1  +  2^7)  ,  (26) 

where  we  have  used  the  subscript  G  for  Gurney.  Since  for  «  =  ®  >  0,  V  >  V_, 

u  G 

Henry  concluded  that  the  jelly  roll  ought  to  give  a  higher  launch  speed  for 
the  outermost  metal  layer  than  the  simple  cylinder  for  the  same  « .  Or,  to. 
achieve  the  same  speed,  V  =  V„,  we  need  n  «  M/C  =  2<tn  =  2M/C~,  that  is  C  = 

u  u  u 

Cq/2,  or  half  the  explosive.  This  kind  of  prediction  is  traceable  to  his 

integrating  over  the  metal  layers  as  if  they  were  expanding  gas,  leading  to 
2  2 

1/4  MV  instead  of  1/2  MV  .  His  use  of  R  instead  of  rg  has  no  practical 

consequences  since  it  cancels  in  Gurney's  scheme  anyway.  If  Henry’s  formula 
were  correct,  it  should  reduce  to  Gurney's  when  all  the  metal  is  placed 
outside  of  the  explosive,  since  he  specifies  nothing  about  ordering  in 
Equation  (92).  Obviously  it  does  not.  Although  Henry's  model  is  incorrect, 
he  has  raised  an  interesting  question  about  the  potential  advantages  of  many 
alternating  layers.  For  example,  instead  of  the  usual  expanding  tube  with 
some  spread  in  wall  thickness  as  air  resistance  slows  the  smaller  fragments 
more  than  the  larger  ones,  we  might  obtain  different  launch  speeds  initially 
which  (together  with  spreading  due  to  air  resistance)  could  fill  more  of  space 
over  a  longer  time  than  a  simple  cylinder  can.  Of  course,  the  slower 
inner-layer  fragments  might  be  less  effective. 

B.  The  Jelly  Roll  and  Onion 

The  first  practical  problem  in  making  such  devices  is  that  of  detonation. 
For  a  head-on  encounter  between  detonation  wave  and  metal  surface,  we  would 
probably  have  to  drill  holea  of  sufficient  size  and  frequency  in  the  metal 
layers  and  fill  them  with  explosive  in  order  to  insure  propagation  to  the 
outermost  layer.  Initiation  in  a  laboratory  device  should  not  be  difficult. 

A  small  spherical  detonator  or  an  exploding  wire  or  foil  might  be  used. 
However,  Inadequate  electrical  sources  might  prevent  the  use  of  such 
techniques  in  field  devices.  Initiation  by  a  disk  or  sheet  of  explosive 
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covering  one  or  both  ends  of  a  jelly  roll  or  Dagwood  should  be  practical  for 
side-on  encounter.  End  initiation  of  a  uniform  explosive  loading  should  lead 
to  an  approximately  plane  wave  consisting  of  the  sum  of  the  waves  in  each 
layer  of  a  jelly-roll  or  Dagwood.  Consequently,  the  pressure  and  density  at 
any  distance  behind  this  front  should  be  approximately  uniform  over  the 
gaseous  portions  of  the  front  and  depend  only  on  the  state  of  expansion.  If 
the  metal  layers  are  significantly  different  in  their  bending  properties,  this 
picture  might  have  to  be  modified.  For  the  present,  let  us  assume  an 
arrangement  for  which  this  picture  is  sufficiently  accurate. 

First,  let  us  consider  the  jelly  roll  and  adopt  the  layer  numbering  scheme 

t  h 

shown  in  Figure  10.  The  specific  mass  (per  unit  length)  for  the  i  layer  of 
explosive  is 

Ci  =  >oiV  (R2i  -  R2i-1  ) 


while  that  for  the  metal  is 

Mi  =  ‘’Mi"  (R2i+1  -  R2i  )  . 

Figure  10  gives  some  examples  for  1=0,  1,  ....  Because  of  cylindrical 
symmetry,  all  R^  are  positive  (R^  =  0  if  1  <  0)  and  are  measured  from  the 

symmetry  axis.  Capital  R^  are  initial  values.  As  the  device  expands,  r^  >  R^ 

are  the  radii,  and  the  decreasing  density  of  each  gas  layer  is 


If  we  assume  that  the  radial  gas  velocity  is  a  simple  linear  form  in  r5  we 
obtain  an  energy  balance  (and  thus  a  launch  speed)  which  is  rather 
complicated.  Thus,  if 


( 91* ) 


(95) 


v  =  A.r  +  B. 

l  l 


\  r2i 


Vir2i-1  \ 
r2i-l  / 


(97) 


t  h 

for  the  i  layer  of  gas,  then  v=V^  when  r=r  and  v  =  V  ^  when  r=r^  1* 

That  is,  the  gas  moves  at  the  metal  speed  V  at  each  interface.  If  1  =  0  is 
the  only  value  of  i,  we  have  v  =  VAr/r)  =  vD(r/rn)  since  r  ,  =  V  ,  =  0,  and 

Equation  (97)  reduces  to  Gurney’s  assumption  for  the  simple  cylinder,  Equation 
(15).  The  energy  balance  becomes 
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c0=  Poo  w  [Ro  " 

C i  "  p0]  7 r  [r2  "  Ri  ] 

^Mo  w  [R]2  "  R°2] 
Ml  ~Pm  !  ^  [^3  '  ^2  ] 


Figure  10.  Layer  Numbering  Scheme  for  the  Jelly  Roll  and  Onion. 
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I  E.C.  =  Js  E  M.V.2  +  H  I  f  Tli  p.v22Trrdr 
1  '  ^i-l 

=  %  L  MV.2  +  2  7  - .  C-~ _  f  21 

1 


2  . 
v  rdr 


If  we  use  Equation  (97)  in  Equation  (98),  we  find 

•  Eici =  H  i  Mivi2  *  ?  i  h  K2  * !  vivi.i  *  <1  ] 

F  =  [i  +  y  x(x2  +  x  -  5)]  /  [l 


+  x(x  -  x  -  1) 


] 


[ 


1  +  x(3x  -  5x 


*!)]/[ 


1  +  x (x^  -  X  -  1) 


] 


(98) 


(99) 

(100) 

(101) 


with  x  =  r^  -\^r21‘  Slnce  x  is  a  ratio  of  radii  at  launch,  we  still  have  no 

dependence  on  device  size.  If  1  =  0  is  the  only  i  value,  Equation  (99) 
reduces  to  Equation  (25). 

If  we  want  to  simplify  our  result,  we  might  follow  Gurney’s  method  of 
choosing  a  form  for  v  designed  for  this  purpose.  Thus,  if 


v 


2 


2 

a.  r 
1 


(102) 


we  still  have  v  =  for  r  =  v  =  ^  for  r  =  r^  ^  and  v  =  VQ(r/ro) 

if  i  =  0  is  the  only  value  of  i.  Since  vdv  *  a^rdr,  the  integral  in  Equation 
(98)  is  simplified  if  we  change  the  variable  to  v  with  limits  1  and  1^. 

By  using  Equation  (102)  in  Equation  (98)  we  find 


E 

i 


E.C. 
1  1 


H  2  M.V. 
xi 
1 


TEci(vi2*Vi-l) 


(103) 


which  also  reduces  to  Equation  (25)  if  i  =  0  is  the  only  value  of  i. 

Now  let  us  consider  the  case  of  a  metal-cored  cylinder  with  C  =  V  =0, 

o  o 

Mq  i  0  and  i  =  0,  1.  Equation  (103)  gives  us 

V1  =  \f(2E1)/(<»1  +  1/2)  .  (104) 

In  other  words,  this  model  predicts  that  a  cored  cylinder  will  have  the  same 
velocity  as  a  simple  cylinder,  if  E1  and  are  the  same.  Of  course,  if  we 


54 


simply  replace  explosive  by  a  metal  core,  we  Increase  and  decrease  V^. 

However,  if  we  alter  the  device  dimensions  in  order  to  keep  a.  the  same,  the 
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launch  velocity  ought  to  be  lower  when  an  inert  core  is  present.  Sterne 

obtained  essentially  the  same  result  by  a  different  method  and  compared  this 
prediction  with  experiment.  For  a  steel  core  with  a  radius  about  3/4  the 
initial  radius  of  an  outer  steel  case,  using  composition  C3  explosive  and  a = 

2.33»  the  measured  launch  velocity  was  down  13Jf  compared  to  that  measured  for 
a  simple  cylinder  with  the  same  and  E^.  He  attributed  this  to  energy 

absorption  and  conversion  to  heat  by  the  immobile  core.  This  hypothesis  was 
approximately  confirmed  by  a  calorimeter  measurement  of  the  heat  content  of 
the  recovered  core  right  after  the  experiment.  He  suggested  lowering  E^  to 

about  Ej  -  10  aQ,  with  E1  expressed  in  calories  per  unit  specific  explosive 

mass.  Here  10  calories  per  gram  of  specific  metal  core  mass  was  approximately 
the  measured  value.  Thus,  for  example,  if  «  700  cal/gm  (as  for  TNT),  then 

a  -  7  would  mean  about  a  10)t  reduction  in  available  energy.  However,  if  «  = 

o  o 

1,  there  would  only  be  a  1.5Sf  reduction.  For  =  1,  =  0.5  in  Figure  8, 

4>  / E  would  be  (0.410)  instead  of  (0.416),  still  larger  than  its  value  (0.333) 

without  a  backing  plate  (aQ  =  0).  Sterne  also  noted  that  there  was  some 

permanent  deformation  of  the  core  which  he  took  to  be  negligible  by  comparison 
to  the  above  correction.  The  energy  absorbed  by  heating  and  fracturing  the 
outer  metal  tube  has  already  been  accounted  for  in  the  original  reduction  of 
E1  from  its  expected  chemical  value  to  its' Gurney  value.  Even  for  a  large 

core  like  that  used,  the  loss  in  velocity  is  fairly  small,  probably  because 
there  is  little  time  for  energy  transfer  by  shock  transmission  and  heat 
conduction. 

As  we  have  noted,  Equation  (99)  is  more  complicated  than  Equation  ( 1 03 ) • 
In  the  case  of  a  cored  cylinder.  Equation  (99)  gives 

V1  =  ^(2E1)/(«1  +  .5F)  ,  (105) 


with  x  =  r^/r^  =  R^/r^  ^  the  radius  of  the  core  does  not  change  during  the 

expansion.  Of  course  for  the  inner  radius  of  the  case  at  launch  r  >  R  >  R„ 

2  2  1 

so  x  <  1.  If  x  =  0  because  =  0  (no  core),  F  =  1  by  Equation  (100)  and 

Equation  (105)  is  that  of  a  simple  cylinder.  Typically  x  might  be  near  0.5. 

As  x  increases  from  zero  to  unity  in  Equation  (100),  F  decreases  from  1  to 
2/3,  being  about  0.8  near  x  =  1/2.  Thus  for  the  same  and  a ,  Equation 

(105)  predicts  that  V1  for  a  cored  cylinder  is  greater  than  for  a  simple 

cylinder,  with  larger  cores  giving  greater  V^.  This  is  contrary  to  the 

experimental  results  discussed  by  Sterne.  Consequently,  Equation  (99)  and 
Equation  (97)  upon  which  it  is  based  are  suspect.  Equation  (103)  and  Equation 
(102)  on  which  it  is  based  are  preferred  in  this  Instance,  not  only  for 
reasons  of  simplicity,  but  also  because  the  prediction  is  at  least  closer  to 
experiment  and  gives  no  misleading  trends. 
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Now  consider  Equation  (103)  for  the  simplest  Jelly  roll  which  is 
illustrated  in  Figure  10  (i  =  0,  1): 

E  C  +  EC  =  1/2  M  V  2  +  1/2  M  V  2  +  1/H  tC  V  2  +  C  (V  2  +  V  2) ]  .  (106) 

o  o  11  oo  11  oo  1o  1 

Equation  (106)  of  course  reduces  to  the  metal-cored  cylinder  if  C  =  V  =0# 

o  o 

Mq  i  0  and  the  simple  cylinder  if  CQ  =  VQ  =  MQ  =  0.  If  only  C1  =  0,  VQ  = 

provided  we  weld  the  two  metal  tubes  together  and  again  we  have  a  simple 

cylinder.  If  the  metal  tubes  are  in  loose  contact  Initially,  the  shock 
rarefaction  wave  returning  from  the  outer  metal/air  interface  will  tend  to 
separate  the  tubes,  while  the  gas  pressure  inside  will  tend  to  keep  them 
together.  Because  of  the  high  gas  pressure  during  the  acceleration  phase  of 
the  motion,  we  expect  the  tubes  to  remain  in  contact  until  they  break  in  most 
cases. 


In  designing  a  warhead.lt  is  desirable  to  have  VQ  <  in  Equation  (106) 

since  this  would  mean  that  the  inner  layer  of  fragments  would  travel  more 
slowly  than  the  outer  layer  and  there  would  be  a  greater  spread  of  metal  in 
space  for  a  longer  time.  This  would  increase  the  chance  of  hitting  a  moving 
target  and  (provided  VQ  and  are  not  too  small)  would  Increase  the 

probability  of  damage. 

It  is  possible  but  not  desirable  to  select  the  E^,  and  in  Equation 

(106)  so  that  the  launch  speeds  are  the  same,  namely,  VQ  *  =  V.  In  this 

case,  with  Eq  =E1  =  E, 


V  =  yj  2E/[ a  +  1/2(1  +  Z> ]  ,  (107) 


where  a  =  M/C  -  (Mq  +  M^)/(Cq  +  C^)  and  Z  =  C^/Q.  If  we  divide  this  equation 
by  the  simple  cylinder  formula. 


V/VQ  =  ^  (a  +  1/2)/[<*  +  1/2(1  +  Z )  ] 


(108) 


for  the  same  a  and  E.  For  Z  >  0  we  have  V  <  VQ  and  neither  an  increase  in 

metal  distribution  is  achieved  nor  an  increase  in  launch  velocity.  For  a 
much  above  unity  V/Vc-+  1.  In  the  extreme  Z  -*  1 ,  (CQ->  0,  VQ-*  0),  which  is 

the  metal-cored  cylinder,  we  can  no  longer  assume  EQ  =  E,  since  the  inner 

metal  tube  absorbs  a  significant  fraction  of  the  energy  Imparted  to  it  and 
converts  it  to  heat  rather  than  to  its  diminishing  kinetic  energy. 
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In  our  comparisons  with  a  simple  cylinder,  we  always  fill  the  simple 
cylinder  with  the  most  energetic  explosive  available  to  avoid  making  the  jelly 
roll  appear  better  because  it  has  a  more  energetic  fill.  We  can,  of  course, 
use  a  less  energetic  explosive  in  either  layer  of  the  simple  jelly  roll. 
However,  the  practical  range  of  E  is  rather  limited.  The  most  energetic  solid 

2 

explosive  commonly  available  is  HMX  with  E  =  (3)  /2,  while  the  least  energetic 

2  _ 

is  TNT  with  E  =  (2.4)  /2.  The  values  in  parentheses  are  the  \/2 E  (mm//ns) 

Gurney  velocities  given  by  Kamlet  and  Finger.  Roughly  speaking  then,  2/3  < 
Eq/E1  <  1  when  EQ  <  E1 .  If  VQ  *  V1  =  V  and  <r  =  EQ/E1  <  1,  then  Equation  (108) 

is  replaced  by  • 

V/Vc  =  V(°  +  [o  +  Z(1  -a)] /[a  +  Ml  +  Z)]  (109) 

so  that  V  will  be  even  less  than  in  Equation  (108),  as  expected.  If  the 
launch  speeds  are  approximately  equal,  then  there  is  not  much  decrease  in 
launch  velocity  compared  to  a  simple  cylinder.  However,  there  is  not  much 
increase  In  metal  distribution  either,  since  the  inner  layer  will  tend  to 
maintain  its  original  separation  from  the  outer  layer  at  least  after  launch 
time.  Most  likely,  the  outer  layer  will  break  first  and  tend  to  slow  down 
first,  emphasizing  the  tendency  for  both  layers  to  fly  together. 

There  are  two  ways  to  achieve  a  significant  increase  in  separation  of  the 
two  layers  compared  to  their  initial  separation.  One  way  is  to  have  MQ/C0  >> 

M^/C^,  approaching  the  case  of  the  metal-cored  cylinder  for  which  Vq  =  0.  The 

other  way  is  to  use  a  less  energetic  explosive  in  the  core  so  that  the  inner 
pressure  is  lower  (at  least  initially).  A  combination  of  both  techniques 
should  enable  us  to  design  a  device  where  the  inner  layer  flies  with  half  the 
speed  of  the  outer  layer,  giving  an  increasing  gap  between  the  layers,  a 
configuration  which  may  be  desirable  for  some  purposes.  An  estimate  of  inner 
layer  launch  speed  might  be  made  by  treating  the  layers  M  ,  CL  and  M.  as  a 
single  Inertial  layer,  so  that  °  1 

Vq  s  ^(2Eo)/[(Ho  +  C1  +  M, )/CQ  +  1/2]  .  (110) 

More  complicated  Jelly  rolls  can  be  discussed  in  the  same  way. 

The  spherical  analog  of  the  jelly  roll  will  also  lead  to  fairly 
complicated  expressions  if  we  use  Equation  (97)  for  v.  However,  if  we  use  the 

cubic  analog  of  Equation  (102)  and  relate  v  to  r  ,  we  find 


O 


I  E.C. 

l  l 


i 


h  Z  M.V. 

ii 

l 


(111) 
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since 

pi  =  MH44  -  4-i ) 

(112) 

(113) 

and 

Mi  -  -  4)  - 

(114) 

using  the 

layer  numbering  scheme  of  Figure  10.  Thus  Equation  (ill)  is 

Independent  of  device  size  and  reasonably  simple.  If  1  =  0  is  the  only  value 
of  i,  we  obtain  Gurney's  sphere  formula.  Equation  (28).  The  same  design  goals 
and  limitations  apply  to  the  onion  as  to  the  jelly  roll,  and  an  entirely 

analogous  discussion  could  be  carried  out.  In  the  limit  of  equal  launch 

o 

speeds,  the  ratio  in  the  last  term  of  Equation  (111)  is  5/3  V  and  Equation 
(ill)  becomes 


EEC.  =  h  (Z  M.  +  L  C.\v 
i  Vi  i  -  h(>4 


N+C)V 


(115) 


For  a  uniform  explosive  fill  of  the  same  type  as  in  a  simple  sphere  with  the 
same  a  =  M/C 

V/Vg  =  4  (a  +  3/5)/(«  +  1)  ,  (116) 

where  V  <  V  ,  the  launch  velocity  of  a  simple  sphere.  A  discussion  of  metal 
s 

spreading  by  controlling  the  E^  and  analogous  to  that  given  for  the 
cylinder  can  also  be  carried  out. 


C.  The  Dagwood 


Here  the  possibility  of  asymmetry  about  the  central  plane  exists, and  the 
momentum  balance  equation  is  not  identically  zero  except  in  the  symmetrical 
limit.  Similar  considerations  hold  for  the  Dagwood  and  jelly  roll  with 
respect  to  detonation  wave  Initiation  and  propagation  as  well  as  pressure  and 
velocity  equalization  by  launch  time.  The  layer  numbering  scheme  In  Figure  11 
enables  us  to  write 


-H 

O 

Q. 

II 

O 

_R2i+l 

R2i]  , 

(117) 

Mi  ■ 

[hi.2  - 

R2i+l]  , 

(118) 

pi  ■  V 

j2i+1  " 

r2i]  . 

(119) 

where  =  Y^  are  initial  values  and  r^  =  y^  are  values  during  the  expansion. 

In  the  planar  case,  using  a  linear  relation  between  v  and  r  enables  us  to 
preserve  Gurney's  basic  assumption.  We  let 
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Cos(/>0)0[Ri-R0]  s  (/>o)o[Ri+lRol] 

Cl=(/>o)l  [R3"  R2] 

C-i  -  {/>o)-i  [R.rRj  s(^o)-i  [|R-2^“^R-t 

Mo  =  (p^)o  [^2"Rl] 

Mis(/>^)i  ^4^3] 

[Ro  **  R-i ]  =  (/>M)-i  [•  R-i  1-lRol] 

Rj  5  0  for  i  <  0  and  Ro  <  0 


Figure  11.  Layer  Numbering  Scheme  for  the  Dagwood. 
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.  v  =  a 


*  b. .  r  1r  *  toa  - 

1  1  L  r2 i-*- 1  '  r2iJ  L  r2i+l  "  r2i  J 


(120) 


since  each  plate  has  a  single  V  at  launch  time  which  may  be  specified  by 
either  the  outer  or  inner  surface  or  some  point  between.  Now  we  find 

Z  E.C .  =  J*  Z  M.V. 2  +  4  z  /vi2  +  v;Vi  i  +  V?  i  ) 

i  i  i  i  1  1  6  i  1  \  1  1  1-1  1-1  '  ( 


(121) 


when  we  use  Equations  (119)  and  (120)  to  find  / 1/2  p,v  dr  and  divide 
3  3  "  1 

V1  -  vi  _  1  by  ai^r2i  +  1  “  r2i^  =  ^vi  "  vi  -  i^*  A  3imllar  Procedure 

involving  Cp^t dr  gives  the  conservation  of  momentum  equation 


v  dr  and  divide 


0  =  Z  M.V.  +  h  Z  C.  (v.  +  V.  «  ) 
i  1  1  i  1 V  1  1-1 ) 


(122) 


which  is  identically  zero  only  for  a  symmetric  Dagwood.  In  general,  Vi  >  0 
for  i  £  0  and  <  0  for  i  <  0.  For  the  simple  asymmetric  sandwich  these 
equations  reduce  to  Equations  (80)  and  (81),  namely, 

E,c,  -  is  [m,v  2  *  m  v  21  +  i  c.  Tv  2  ♦  v.v  *  v  2 1 
11  Ln  oo  J  6  1  Ll  1  o  o  J  (1. 

and 

0  =  Tm.v.  +  M  V  1  +  h  c.  Tv.  +  V  "1 

L  1  1  °  °J  1  L  1  °J  (12 


(123) 


-  |_“r  1  o  oj  --1  L'l  Oj  (12U) 

if  we  identify  E1  =  E  and  C1  s  C  (with  CQ  *  M  1  =  C_,|  =  M_2  =  0  in  Figure  11). 

Since  for  given  M^,  and  E  we  have  two  equations  in  two  unknowns,  we  can 

solve  for  V  and  V.  as  before, 
o  1 

For  the  asymmetric  Dagwood  illustrated  in  Figure  11,  Equations  (121)  and 
(122)  become 


E-1C-1  +  EoCo  +  E 


ici  -  ■»  [M. 


,V  ,2  +  M  ,V  ,2  +  M  V  2  +  M 
2-2  -1-1  o  o 


iV] 


i|C-l  (V-l2  *  V.lV-2  ♦  V-22) 


(125) 


+  c  (v2  +  vv.+v2) 

o  \  o  o-l  -1  / 

*  C1  (V l2  *  Vo  *  Vo2  )  I 

0  =  [m_2V_2  *  M_JV_1  *  McVo  ♦  MjV,] 

^{C-l(V-l*V-2)*Co  (Vc  +  V-l)tCl(Vl*Vo)} 


(126) 
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Now  for  given  and  we  have  two  equations  In  four  unknown  V.^.  For 

the  symmetric  case  when  Equation  (126)  Is  identically  zero,  we  have  one 

equation  In  two  unknowns,  v  =  -V,  and  V,  =  -V«. 

o112 

Let  us  consider  the  symmetric  case,  assuming  equality  of  velocities  at 
launch  time  and  a  single  explosive  loading.  With  er  =  2(MQ  +  M1)/(Cq  +  C^)  = 

M/C,  we  find  that  Equation  025)  with  E^  =  EQ  =  E^  =  E  becomes 

f 

2E  s  («+  1)V2  (127) 

for  the  simplest  symmetric  Dagwood  with  explosive  in  the  center.  This  may  be 
compared  with  Equation  (83)  and  V  for  the  simple  symmetric  planar  case  to 
obtain  (for  the  same  E  and  «)  P 

V/Vp  s  y]  (a  +  1/3)/(«  +1)  028) 

so  that  the  symmetric  Dagwood  launch  speed  Is  somewhat  less  than  for  the 
simple  symmetric  sandwich.  For  typical  a  values  the  loss  in  launch  speed  will 
not  be  great,  but  for  velocity  equalization  the  Increase  in  metal  spreading 
will  not  be  great  either.  For  the  ideal  planar  configuration  (unlike  the 
cylinder)  there  will  be  no  metal  layei*  thinning  or  fragmentation  during 

flight.  A.  launch  In  air  instead  of  vacuum  will,  however,  tend  to  slow  the  outer 

plates  first.  Again  we  may  increase  the  spread  of  metal  at  the  sacrifice  of 

inner  layer  speed  by  controlling  the  ratios  E^/E-p  M^/C^  and  M^/C^.  In 

addition,  an  estimate  of  VQ  can  be  made  by  treating  the  layers  exterior  to  CQ 

as  a  single  inertial  mass. 

If  M_2  s  C-1  =  v0  =  M_i  =  Mt»  c0  =  «  =  (2M1)/2C1)  and  = 

V2  In  Figure  11,  we  have  a  symmetric  sandwich  with  the  metal  plate  Mq  at  Its 
center.  For  EQ  =  =  E,  Equation  (125)  becomes 

V  =  ^(2E)/(«  +  1/3)  ,  (129) 

and  the  model  predicts  the  same  launch  speed  as  If  M  =  0  (the  simple 

symmetric  sandwich).  Again  we  should  account  for  the  fact  that  the  effective 
values  Eq  s  E^  <  E  since  some  energy  Is  absorbed  by  the  Inert  stationary  core 

as  In  the  case  of  the  metal-cored  cylinder. 

If  C  j  s  C  s  0  in  Figure  11  and  Equations  (125)  and  (126),  we  have  four 

plates  but  only  one  layer  of  explosive  in  the  middle.  Again  we  have  a  simple 
asymmetric  sandwich,  provided  we  bond  the  plates  together  so  they  will  not 
trap  momentum  and  separate  as  soon  as  rarefaction  waves  reflected  from  the 
free  surfaces  return  to  the  Interfaces.  This  configuration  is  a  generaliza¬ 
tion  of  the  case  discussed  by  G.  Jones'®  which  consists  of  more  than  one 
plate,  but  only  on  one  side  of  a  single  sheet  of  explosive. 
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If  asymmetries  are  not  great,  and  we  assume  approximate  equalization  of 
velocities  at  launch.  Equation  (127)  holds  for  any  number  of  explosive  sheets 
sandwiched  between  metal  plates. 

For  significant  asymmetries,  we  encounter  a  problem  as  soon  as  we  add  a 
sheet  of  explosive  to  one  side  of  a  simple  asymmetric  sandwich.  If  we  let  E 

o 

=  E1  s  E,  H  2  :  H  1  =  C  1  =  0  in  Figure  11,  Equations  (125)  and  (126)  become 

EC  =  1/2[M  V  2  +  M.V  2]  .  1/6[C  (V  2  +  V  V  .  +  V2  ,)  .  C,(V  2  +  V,V  *  V  2) ] 

0  0  11  o  o  o-l  -1  11  10  o(130) 

and 

0  =  [MqV0  +  M^]  +  1/2[Cq(V0  +  V_,)  +  C1(V1  +  V0)]  ,  (13D 

where  V  1  is  the  launch  velocity  of  the  free  gas  surface,  which  cannot  be  zero 

if  CQ  >  0  in  the  present  configuration  (M^/C1/MQ/C0) .  Since  we  only  have  two 

equations  for  the  three  unknowns  V.,  VQ  and  V  , ,  we  cannot  give  a  general 

solution  within  the  limitations  of  this  model.  If  for  some  reason  we  want  to 
fix  the  V^,  then  we  would  have  two  linear  equations  in  some  M^,  pair  (for 

given  E)  and  some  physically  interesting  solutions  would  exist. 

The  assumptions  of  uniform  density  and  linear  spatial  distribution  of  gas 
velocity  at  any  instant  imply  the  existence  of  a  stationary  plane  (or  surface) 
of  zero  gas  velocity.  This  was  mentioned  above  in  the  discussion  after 
Equation  (90).  For  the  simple  asymmetric  sandwich  (M1/C1/MQ)  the  location  of 

this  plane  may  be  found  by  letting  v  =  0,  r  =  R,  and  1  =  1  in  Equation  (120), 
giving  R  =  (V^  -  v0R3^/^vi-V0)  with  VQ  =  _qVi  and  Vl  Siven  by  Equation  (82). 

Thus  we  can  find  R  from  a^,  aQ,  R2,  and  R^  since cancels.  For  example, 

when  ao  =  £r.j,  q  =  1,  and  R  =  (R^  +  R^)/2.  That  is,  the  plane  lies  in  the 

middle  of  the  layer  If  we  add  the  layer  CQ  below  the  sandwich  (M1/C1/Mq), 

we  shift  the  plane  downward.  Sufficiently  large  CQ  will  move  the  plane 

through  M  and  into  C  .  When  v  =  0,  we  have 
o  o  o 

V1  sy2E  /  I  ttl  +  1/3  +  py  +  2ay)2n  1  -  j  (132) 

and 

V_1  =“  1  (^i  +  2«1)/d  -/?,)]  vi  ,  (133) 

where  =  C^/C  =  C^/(CQ  +C^)s1-0s1-  Cq/C.  This  case  is  equivalent 

to  the  simple  open-face  sandwich  (M^C)  with  a  reduction  in  E  dependent  on  the 

size  of  Mq  as  explained  above  for  other  "cored"  configurations.  Usually  this 

correction  is  negligible.  If  we  neglect  it  and  set  Equation  (132)  equal  to 
Equation  (84),  we  find  /?1  =  1/[2(1  +  a^)].  Obviously  for  =  0,  p^  =  pQ  = 

1/2.  For  =1,  /3.J  =  1/4  and  Mq  must  divide  C  so  =  C/4  and  Cq  =  3C/4  to 
make  V0  =  0,  etc. 
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If  we  wish  to  launch  two  plates  in  the  same  direction  but  with  V  <  V.  so 

o—1 

there  will  be  a  greater  spread  of  metal,  our  model  is  too  simple  to  enable  us 
to  analyze  all  possibilities.  However,  for  certain  fixed  relations,  we  can 
estimate  what  is  needed.  For  example,  to  obtain  vq  =  sV^  with  s  <  1,  using  MQ 

=  kM^  with  k  >  1,  Equations  (131)  and  (130)  become 

V_1  =  -CE  (s  +  (i^)  +  2ay  ( 1  +  ks)  ]/( 1  -  y31)]V1  (134) 

and 


_  (135) 

For  example,  if  s  =  1/2,  k  =  8,  s  1/2  and  (3^  =  1/4,  we  find  V ^/y/2E  =  0.55 
from  Equation  (135)  with  VQ/\/~2E  =  0.275.  Of  course  we  can  also  use  different 
types  of  explosive  in  the  two  sheets  (Eq  i  E^)  and  vary  the  design  further. 

The  principles  are  clear  and  will  not  be  elaborated  here. 

The  next,  more  complicated  case  is  (M1 /C1 /M0/CQ/M_1 )  in  Figure  11.  This 
requires  the  addition  of  1/2M  ^  ^  to  the  right  side  of  Equation  (130)  and  of 
M  V  to  the  right  side  of  Equation  (131).  There  is  one  more  dimensionless 
parameter,  a  1  =  M  ^/ C ,  but  still  only  three  launch  speeds.  Let  us  close  this 

section  by  giving  some  simple  examples  of'  how  we  might  estimate  launch  speeds. 
Suppose  we  have  =  MQ  =  M^,  CQ  =  2C^,  and  EQ  =  E^.  We  can  solve  the 

momentum  equation  for 

V-1  *  ao  vo  *  a1  V * * * * X1  .  <136) 

where 

a  =  -(3+60/(2+60  (137) 

O  I  If 

and 


a1  =  -( 1+6a1 ) /(2+6a1 )  (138) 

and  use  Equation  (136)  to  eliminate  If  1  in  the  energy  equation  to  obtain 

2E  =  XV  2  +  YV  v  +  ZV  2  (139) 

o  o  1  1  » 


where 


X  =  ac2(«  *2/9)  +  2/9aQ  +  («1+1/3)  (140) 

Y  =  2aQa1(a1+2/9)  +  2/9aQ  +1/9  (141) 

Z  =  a^2(«1+2/9)  +  («1+1/9)  >  (142) 
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Thus  the  V  depend  only  on  a  =  M  /C,  where  C  =  C  +C  .  Suppose  we  let  a  =  1, 

1  11  p  0  1  p 

that  Is,  a  s  M/C  =  31^/C  =  3*  If  2E  =  (2.4)  =  5.76  (mm/^s)  as  for  TNT  at  a 

typical  loading  density,  then  Equation  (139)  becomes 

5.76  =  2.63  VQ2  +  2.27  VQV1  +  2.05  ,  (143) 

while  Equation  (136)  becomes 

V  ,  =  -(1.125V  +  .875V,).  (144) 

-1  o  1 

If  we  consider  two  related  cases  which  are  known,  we  can  restrict  the  range  of 
allowed  values.  First,  consider  Cq  =  Cj  =  C/2,  a  case  for  which  Vq  =  0  and  V1 

=  2.H//2M.]/C  +  1/3  =  1.57  mm/^s  =  -V^*  If  we  move  Mq  so  C^C/2  and  Cq>C/2, 
then  Vq>0,  V^l.57  mm//is  and  V_1<-1.57  mm//is.  Second,  consider  Cq  =  ^  =  C/3, 
for  which  VQ  =  0  and  V1  =  2.4//2M1/(2C/3)  +  1/3  =  1.315  mm//iS  =  -V^. 
Increasing  Cq  to  2C/3  =  2C1  will  make  V^>1.315  mm/fis  and  VQ>0.  The  range  of 
V.|  values  is  thus  fairly  narrow,  namely,  1-31 5< V ^ <  1 . 57  mm//j.s).  These  extreme 
values  may  be  used  in  Equation  ( 1  il 3 )  to  find  the  range  of  VQ,  namely 
.51 2>Vo>.  176  and  from  Equation  044)  we  find  -1 .727>V_1>-1 .572.  If  we  take 

the  mid-points  of  these  ranges  as  our  estimate,  we  see  that  the  model  predicts 
V i  =  1.44+.13  mm//is,  VQ  =  .34  +.17  mra//is  and  V  ^  =  —  1 . 6i< . 08  mm/MS. 

Suppose  we  wish  to  make  V1  =  2Vq>0.  Of  course  there  are  many  choices  of 

M.  and  C,  which  will  accomplish  this  goal.  For  example,  a  -  «,  =  1/2,  «  ,  = 
11  oi-i 

1  and  a  s  M/C  =  (a,+<*  +a  ,)  =  2.  In  this  case.V,  =  1.5*1  mm/^s,  V  =  .77  mm//*s 
1  o  -1  1,0 

and  V  =-1.76  mm//is  for  =  2.4  mm/MS.  In  this  and  the  previous  example  we 
are  neglecting  any  reduction  in 

The  procedure  illustrated  here  is  not  completely  straightforward  since  it 
requires  some  ingenuity  to  narrow  the  range  of  values  by  comparisons  with 
known  theoretical  or  experimental  cases.  However, it  can  give  us  useful 
estimates  for  fairly  complicated  Dagwoods,  jelly  rolls  and  onions. 
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SUMMARY 


In  this  report  we  have  seen  that  the  exponent  V  In  the  entroplc  equation 
of  state  may  be  taken  equal  to  three  in  practical  applications  involving  the 
expansion  of  solid  explosive  detonation  products.  This  fact  has  been  used  to 
simplify,  extend  and  unify  the  ideas  of  Taylor,  Gurney,  Thomas  and  Sterne 
concerning  the  acceleration  of  metal  by  explosives.  The  modified  theories 
which  have  been  developed  here  can  be  used  in  the  design  of  lined  cavity 
charges  and  fragmentation  warheads  as  well  as  flying  plate  experiments  like 
those  which  have  been  used  to  study  the  shock  properties  of  solids. 
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director  of  the  Ballistics  Research  Laboratory  who  advocate  a  return  to 
fundamentals  as  a  means  of  obtaining  fresh  Insights  and  of  stimulating  new 
approaches  to  problems.  The  author  shares  this  view  and  has  attempted  to  give 
it  substance  in  this  report. 
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